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The  objectives  of  this  program  were  to  develop  a  new  growth  technique  for  _ 
indium-doped  silicon,  to  determine  the  trapping  parameters  for  the  centers 
in  this  material,  and  to  identify  and  if  possible  eliminate  the  X-center,  j 
an  indium-related  complex.  j 
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The  growth  technique,  studied  under  this  contract  was  the  solution  growth  i 
of  silicon  from  an  indium  melt.  Crystals  doped  with  1.6  x  1013  In/cm3  i 
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were  grown  by  this  technique  at  'v-1350°C.  This  is  at  the  maximum  in  the 
Indium  concentration  versus  temperature  curve  and  represents  a  concentra¬ 
tion  three  to  five  times  the  values  obtained  from  Czochralski  or  float-zone 
growth.  The  wafers  grown  have  a  uniformity  of  better  than  +  5%  by  spread¬ 
ing  resistance  both  axially  and  radially.  Boron  and  phosphorus  can  be  held 
to  below  5  x  10 12  atoms/cm'  with  the  boron  compensated  to  less  than  108 
electrically  active  boron  atoms/cm  .  The  indium- X-center  concentrations 
can  be  held  below  the  detection  limit  of  about  5  x  1013/cm.  One-inch- 
diameter  wafers  with  no  voids  or  polycrystalline  areas  were  grown  but  at 
present  larger-diameter  growths  suffer  from  these  defects.  The  current 
solution-growth  dislocation  count  is  'v-1000  counts/cm2. 

The  trapping  parameters  of  indium  and  the  indium-X  center  were  measured  as 
a  function  of  temperature  and  electric  field  by  deep-level  transient  spec¬ 
troscopy.  Both  centers  show  the  strong  field  dependence  expected  for 
coulombic  traps.  Two  deep  hole  traps  at  E  +  0.35  eV  and  E  +0.41  eV  were 
found  in  indium-doped  silicon. 

X-centers  were  found  associated  with  all  of  the  other  group  IIIA  dopants  by 
optical  absorption  and  deep-level  transient  spectroscopy  techniques.  The 
defects  depend  on  carbon,  and  their  annealing  kinetics  fit  the  reaction 
acceptor  +  carbon  j  X.  The  aluminum-X  center  has  the  stress  symmetry  of  a 
|lll[  trigonal  defect.  It  was  concluded  that  the  best  model  for  these 

defects  is  a  nearest-neighbor  carbon-acceptor  pair.  By  growing  with  low- 
carbon  materials,  the  X-center  can  be  eliminated. 
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SECTION  1 
INTRODUCTION 


Areas  of  current  interest  in  infrared  detector  development  include 
the  development  of  larger  arrays,  the  use  of  cheaper  materials, 
and  the  integration  of  the  detectors  and  processing  electronics 
onto  the  same  chip.  Silicon  doped  with  various  impurities  can  be 
used  as  a  photodetector  in  many  regions  of  the  infrared.  Silicon 
is  available  in  large  areas  and  the  technology  and  equipment  need¬ 
ed  to  process  the  detectors  and  coupled  electronics  is  that  used 
in  the  current  silicon  electronics  industry.  Because  of  this  the 
production  of  chips  with  integrated  detectors  and  electronics 
would  be  cheaper  with  silicon  than  with  any  other  material.  The 
yields,  uniformity,  and  quality  of  silicon  devices  might  be  ex¬ 
pected  to  be  better  for  silicon  than  for  any  other  less  highly 
developed  material  and  the  development  time  and  cost  should  be 
much  less. 

This  report  deals  with  the  development  of  indium-doped  silicon  as 
a  detector  material. 

Detector-grade  indium-doped  silicon  should  have  high  indium  con¬ 
centration,  good  uniformity  over  2-inch-or-larger  wafers,  and  a 
minimum  concentration  of  electrically  active  shallow  defects.  A 
shallow  defect  complex  which  has  been  found  in  indium-doped 
silicon  is  the  indium-X  center.  This  center  along  with  uncompen¬ 
sated  boron  can  introduce  thermal  noise  into  the  detector  and 
require  a  lower  temperature  of  operation  than  is  needed  if  these 
defects  are  absent.  In  some  applications  the  extra  cooling  re¬ 
quired  when  X-levels  are  present  results  in  a  severe  handicap. 
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The  specific  aims  of  this  study  effort  were  three  fold: 

1)  Develop  a  new  growth  process  which  can  produce  better- 
quality  indium-doped  silicon  than  current  Czochralski 
or  float-zone  techniques. 

2)  Determine  the  origin  of  the  X-center  and  if  possible 
eliminate  it. 

3)  Determine  the  trapping  parameters  of  centers  in  indium- 
doped  silicon 


This  report  reviews  the  results  of  the  total  study  effort  with 
emphasis  on  the  results  of  the  last  reporting  period.  Many  of  the 
results  are  now  published.  Results  which  are  not  yet  in  the  open 
literature  are  included  in  the  appendices. 
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SECTION  2 
SOLUTION  GROWTH 


DESCRIPTION  OF  THE  TECHNIQUE 

Solution  growth  is  a  technique  for  growing  silicon  crystals  from 
a  solution  of  silicon  in  molten  indium.  The  apparatus  is  shown 
in  Figure  1.  The  materials  used  are  0.999999  pure  indium  and  the 
vacuum  float-zone  20,000  ft-cm  silicon.  A  sandwich  is  formed  with 
a  silicon  source  crystal  on  top,  an  indium  layer,  and  a  seed 
crystal  on  bottom.  The  thin  indium  layer  and  the  wafer-sized 
source  and  seed  crystals  provide  for  a  minimum  of  material  wasted 
in  the  growth  process.  At  high  temperatures  silicon  dissolves  in 
the  indium.  A  temperature  gradient  with  the  source  hotter  than 
the  seed  produces  a  silicon  concentration  gradient  in  the  solution, 
with  more  silicon  dissolved  near  the  hotter  source  crystal .  This 
concentration  gradient  moves  silicon  from  the  source  to  the  seed 
where  it  recrystallizes.  The  recrystallized  material  contains  the 
maximum  indium  concentration  allowed  at  the  growth  temperature. 

The  amount  of  material  grown  depends  on  the  growth  temperature  and 
the  temperature  gradient.  At  1300°C  'v  10  mm  of  silicon  can  be 
grown  in  24  hours  with  a  1.7°/mm  temperature  gradient.'*'  The  growth 
rate  and  the  physics  of  the  solution  growth  technique  are  dis¬ 
cussed  in  detail  in  Appendix  A. 


ht.  Scott  and  R.J.  Hager,  J.  of  Electron.  Mater.  8,  581  (1979) 
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Figure  1.  Schematic  of  the  close-spaced 
solution-growth  apparatus. 


DOPANTS  IN  SOLUTION-GROWN  MATERIAL 

The  solubilities  of  indium,  oxygen,  and  carbon  as  a  function  of 
temperature  in  silicon  have  been  measured  with  a  series  of  solution 
growths.  The  indium  concentrations  were  determined  by  Hall,  resis¬ 
tivity,  and  optical  absorption  measurements  (see  Appendix  A).  The 
oxygen  and  carbon  concentrations  are  measured  by  optical  absorption 
(see  Appendix  B).  The  data  are  shown  in  Figure  2.  The  solubility 
of  indium  is  retrograde,  with  the  maximum  being  1.6  x  1018  In/cm3 
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Figure  2.  The  maximum  indium,  carbon,  and  oxygen 
concentrations  as  a  function  of  growth 
temperature  as  determined  from  solution 
growths.  The  solid  lines  represent  the 
oxygen  and  carbon  solubility  curves 
obtained  by  Hrostowski  and  Kaiser  and 
Bean  and  Newman.  Data  from  Czochralski 
crystals  are  shown  at  the  silicon 
melting  point  of  1420°C. 


^H.J.  Hrostowski  and  R.H.  Kaiser,  J.  Phys.  Chem.  Solids  9,  214 
(1958). 

^A.R.  Bean  and  R.C.  Newman,  J.  Phys.  Chem.  Solids  32,  1211  (1971). 
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occurring  at  1350°C.  This  concentration  is  about  three  times  the 
best  Czochralski  values  and  five  times  the  float-zone  values  ob¬ 
tained  from  growths  at  the  silicon  melting  point  of  1420°C.  Oxy¬ 
gen  seems  to  be  introduced  from  the  indium,  since  it  is  seen  even 
in  growths  using  low-oxygen  float-zone  silicon  which  has  been  iso¬ 
lated  from  the  quartz  walls.  The  oxygen  concentrations  seen  in 
solution  growths  are  always  approximately  the  solubility  limit  for 
the  growth  temperature.  This  is  seen  in  Figure  2  by  comparing  the 

solution  growth  data  points  and  the  oxygen  solubility  curve  ob- 

2 

tamed  by  Hrostowski  and  Kaiser.  Carbon  has  been  introduced  from 

the  silicon  source  and  from  graphite  added  to  some  growths.  The 

carbon  data  shown  are  compared  with  the  curve  obtained  by  Bean 
3 

and  Newman.  The  maximum  substitutional  carbon  concentration  seen 

17  *5 

in  the  solution  growths  is  in  the  mid-10  atoms/cm  range,  which 
agrees  with  the  more  recent  work  of  Newman3  but  is  a  factor  of 
three  lower  than  that  reported  in  Newman's  earlier  work.^  The  car¬ 
bon  content  can  be  controlled,  and  the  use  of  low-carbon  source 
and  seed  silicon  materials  has  yielded  crystals  with  carbon  con¬ 
centrations  below  our  optical  detection  limit  of  ^5  x  lO'*'3  atoms/ 
cm3 . 

Boron  and  phosphorous  are  introduced  from  the  source  and  seed 
materials.  Boron  can  also  be  introduced  from  the  quartz  holders. 
Solution  growths  made  from  vacuum  float-zone  silicon  typically 
have  boron  and  phosphorous  concentrations  less  than  5  x  1012/cm3 
as  measured  by  optical  absorption  using  a  bandgap  bias  light  to 

populate  the  compensated  boron  and  the  phosphorous  donors.  These 
values  are  very  close  to  those  obtained  from  the  original  source 


4 

R.C.  Newman  and  J.B.  Willis,  J.  Phys.  Chem.  Solids  26,  373 
(1964) 
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and  seed  crystals,  showing  that  the  solution  growth  technique  does 

not  introduce  any  appreciable  concentrations  of  these  impurities. 

Hall  data  show  that  the  boron  concentrations  are  compensated  at 

8  3 

least  down  to  the  10  carriers/cm  by  the  residual  phosphorous. 
These  data  are  shown  in  Figure  3. 


Figure  3.  Hall  data  for  solution-grown  silicon 
showing  indium  dominating  the  free- 
carrier  generation  down  to  about 
the  108/cm3  level. 
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Optical  absorption  was  also  used  to  monitor  the  electrically  active 
aluminum,  gallium,  and  thallium  concentration  in  the  solution- 
grown  material .  The  concentrations  of  all  of  these  elements  were 

•j  2  o 

below  the  optical  detection  limit  of  x  10  /cm  A1  or  Ga  or 
'V10-LO  Tl/crn  . 

Spreading  resistance  measurements  have  shown  the  solution-growth 
doping  to  be  uniform  to  better  than  ±5%  both  axially  and  radially. 
This  is  shown  in  Figure  4,  where  data  from  Czochralski  and  float- 
zone  crystals  are  shown  for  comparison.  The  Czochralski  and  float- 
zone  samples  show  concentration  drift  from  the  center  to  the  edge 
of  a  wafer  and  they  also  show  striations  due  to  the  crystal  rota¬ 
tion  during  growth.  Solution  growth  is  an  equilibrium  process 
carried  out  at  a  constant  temperature  with  no  motion  to  the  ma¬ 
terials.  Growths  near  1350°C  are  also  at  the  broad  top  of  the 
indium  solubility  curve  so  that  small  temperature  variations  make 
little  difference  in  the  indium  solubility. 

The  X-center  in  indium-doped  silicon  has  been  found  to  depend  on 

carbon  and  annealing  temperatures.  ’  The  data  are  discussed  in 

the  section  on  X-center  studies.  Solution  growths  made  from  low- 

carbon,  20,000  fi-cm  vacuum  float-zone  silicon  and  annealed  at 

600°C  for  6  hours  have  shown  no  detectable  X-center  concentration 

down  to  the  optical  detection  limit  of  ^5  x  1013/cm3.  The  Hall 

data  of  Figure  3  show  no  electrically  active  X-centers  down  to  the 
8  3 

10  /cm  level.  Thus,  one  result  of  this  study  has  been  the  elimina¬ 
tion  of  the  indium-X  center. 

5R.  Baron,  J.P.  Baukus,  S.D.  Allen,  T.C.  McGill,  M.H.  Young, 

H.  Kimura,  H.V.  Winston,  and  O.J.  Marsh,  Appl .  Phys .  Lett.  34, 

257  (1979).  — 

ry 

dC.E.  Jones,  D.E.  Schafer,  W.  Scott,  and  R.  Hager,  Bull.  Am. 

Phys.  Soc.  24,  276  (1979). 


8 


QJ 

> 

(— 

3 

LlI 

cc 


UJ 

o 

< 


co 

l/) 

LU 

a: 

o 


Q 

< 

UJ 

QJ 

CL 

CO 


100 


10 


SOLUTION 


C20CHRALSKI 


CENTER 


FLOAT  ZONE 


EDGE 


DISTANCE  (mm) 

Figure  4.  Spreading  resistance  measurements  from  edge  to  center 
of  a  solution  growth,  a  Czochralski  growth,  and  a 
float-zone  growth  for  indium-doped  silicon. 


VOLUME  PRODUCTION 

Solution-grown  indium-doped  silicon  has  shown  promise  because  of 
its  high  indium  concentrations,  good  uniformity,  and  low  impurity 
levels.  Solution-grown  silicon  will  not  be  practical  for  large- 
scale  detector  production  unless  it  can  be  made  in  the  volume 
needed ,  however . 
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A  small  furnace  configuration  has  been  designed  and  tested.  This 
is  shown  in  Figure  5.  The  system  involves  a  cylindrical  furnace 
area  ^5  inches  in  diameter  by  ^2  inches  high  heated  by  a  plate 
heater  on  the  top.  In  the  test  furnace,  a  gradient  of  ^1.5°C/mm 
over  the  temperature  range  of  1000  to  1420°C  can  be  established 
vertically  down  from  the  plate  heater.  This  could  be  maintained 
to  better  than  ±0.1°C.  The  growth  rate  estimated  for  this  gradient 
is  M  mm/ day  at  1200°C  and  10  mm/ day  at  1300°C  (see  Appendix  A). 
Test  growths  have  been  made  between  1100°  and  1200°C  which  fit 
this  prediction.  The  furnace  can  produce  wafers  2  to  3  inches  in 
diameter . 


5" 


Figure  5.  Furnace  and  sample  arrangement  for 
close-spaced  solution  growth. 


The  simple  nature  of  the  furnace  design  and  its  small  size  allow 
a  large  number  of  furnaces  to  be  set  up  in  parallel  without  undue 
cost.  Parallel  growth-furnaces  coupled  with  the  growth  rates  mea¬ 
sured  could  produce  adequate  volumes  of  solution-grown  crystals 
for  detector  needs  . 
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CRYSTAL  QUALITY 


Most  of  the  data  on  solubilities  and  purities  have  been  taken  on 
small  samples  on  the  order  of  2  cm  in  diameter.  To  be  useful  for 
large-array  production,  crystals  2  inches  in  diameter  or  larger 
are  needed.  This  is  the  minimum  size  which  can  be  handled  with 
ease  on  present  silicon-device  production  lines.  The  wafers  should 
be  uniform  in  doping,  single-crystal  and  void-free,  and  have  low 
dislocation  counts.  The  major  effort  in  the  last  period  of  this 
contract  was  directed  toward  scaling-up  the  solution-growth  diam¬ 
eter  while  maintaining  crystal  quality.  Figure  6  is  a  photograph 
which  shows  the  growth  sizes.  The  larger  the  crystal  area,  the 
harder  it  has  been  to  grow  wafers  without  some  voids  or  poly¬ 
crystalline  areas.  Progress  has  been  made  from  1-cm-diameter 
wafers  without  these  defects.  At  present,  growths  larger  than  1 
inch  in  diameter  have  all  had  some  voids  or  polycrystalline  areas. 
A  series  of  solution  growths  have  been  made  to  try  to  determine 
the  parameters  affecting  the  start  of  polycrystalline  growth  and 
void  formations. 

Polycrystalline  areas  can  be  nucleated  from  defects  in  the  seed 
layer,  from  particles  in  the  melt,  and  from  growth  flaws.  The 
quality  of  solution-growth  material  has  been  noted  to  drop  at  the 
height  of  the  original  indium  solution,  presumably  due  to  impur¬ 
ities  which  float  on  the  top  of  the  original  melt.  The  crystal 
quality  for  solution-growth  runs  with  various  indium  melt  depths 

running  from  0.125  inch  to  0.002  inch  was  checked.  Again,  a  ten¬ 
dency  for  polycrystalline  growth  or  complete  crystal  reorientation 
was  found  at  the  position  of  the  top  of  the  original  indium  layer. 
Better  growths  were  obtained  with  thinner  indium  layers,  with  the 
optimum  indium  thickness  being  v0.005  inch. 
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Figure  6.  Progress  in  solution-growth  wafer  size  from  1  cm 
in  diameter  in  1977,  to  1  inch  in  diameter  grown 
in  1978,  to  a  1.75-inch-diameter  wafer  grown  in 
1979. 

Crystal  cracking  has  been  found  when  the  indium  spills  out  of  the 
well  and  fuses  the  bottom  of  the  seed  crystal  to  the  quartz  ampoule. 
The  differential  contraction  on  cooldown  causes  the  cracking.  This 
has  been  prevented  by  placing  the  seed  crystal  on  a  silicon  or 
quartz  ring  support  which  is  smaller  in  diameter  than  the  seed 
crystal . 

Voids  have  been  postulated  to  be  due  to  trapped  gas  in  the  indium, 
indium  vaporization,  indium  spillage,  or  incomplete  indium-silicon 
wetting.  Test  runs  have  been  made  on  pre-wetting,  on  preheating  to 
drive  off  trapped  gases,  and  on  indium  overpressures  to  prevent 
vaporization.  Molten  indium  has  a  strong  tendency  to  ball  up  and 
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not  wet  a  silicon  surface  unless  the  temperature  is  raised  to 
'I'llOO  C.  Solution  growths  are  now  put  in  a  radiant  heater  and 
brought  to  1100°C  for  10  minutes  in  argon  or  vacuum.  The  source- 
indium-seed  sandwich  is  then  X-rayed  to  see  if  any  voids  are 
present.  Even  after  complete  wetting,  voids  have  been  found  to 
form  during  growth.  This  is  especially  true  in  vacuum  growth. 
Indium  vaporization  seems  to  be  a  major  cause  for  void  formation. 
Indium  and  inert  gas  overpressures  have  helped  to  reduce  this 
problem.  Solution  growths  currently  are  done  in  a  small  quartz 
container  which  helps  maintain  the  indium  vapor  in  the  growth 
region.  This  is  shown  in  Figure  5.  Some  voids  are  still  seen  in 
the  larger-area  growths.  Growths  in  saturated  indium  vapor  or 
high-pressured  inert  gas  have  not  been  attempted  yet  but  they  may 
help  to  eliminate  voids. 

Dislocation  counts  have  been  made  for  solution-growth  crystals 

over  various  orientations  of  substrates.  The  seed  material  has 

2 

averaged  10  dislocations/cm  .  Some  growths  measured  within  a  few 
mils  of  the  seed  have  been  as  good  as  the  seed  quality.  Most  solu¬ 
tion-growth  material,  however,  has  shown  dislocation  counts  of 
2 

^1000/cm  independent  of  the  substrate  orientation.  Figure  7  is 
a  photograph  of  an  etched  solution-growth  silicon  surface. 

For  a  more  detailed  discussion  of  solution  growth,  the  reader  is 
referred  to  Appendix  A.  The  material  in  Appendix  A  was  presented 
at  the  June  1979  IRIS  Specialty  Meeting  on  Infrared  Detectors  and 
covers  the  growth  rates  and  impurity  analysis  in  greater  detail 
than  is  found  in  this  section. 
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Figure  7.  Solution-growth  silicon  etched  with  Wright 
etch  to  show  dislocation  and  defects.  The 
etch  pit  count  shown  is  'vlOOO/cm2. 
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SECTION  3 
X-CENTER  STUDIES 


The  X-center  is  an  indium-related  complex  which  is  seen  in  Hall 

7  8 

and  optical  absorption  data  at  0.12  eV.  ’  This  center  is  a  source 
of  noise  in  indium-doped  silicon  detectors  unless  the  device  is 
operated  at  a  temperature  of  40  to  50K  to  freeze  it  out .  The  normal 
operating  temperature  for  an  indium-doped  detector  should  be  60 
to  70K.  The  extra  cooling  needed  when  X-centers  are  present  is  a 
severe  handicap  in  some  applications. 

There  are  several  models  which  have  been  proposed  for  the  X-center. 

5 

Both  the  Hughes  group  of  Baron  et  al.  and  this  group  have  seen 

a  carbon  dependence  for  the  X-center  concentration  and  have  sug- 

6 

gested  a  carbon-acceptor  nearest-neighbor  pair  model.  The  X-center 
concentration  is  seen  to  increase  after  irradiation  and  the  AFML 
group  has  suggested  that  the  center  may  involve  interstitials  or 

9 

vacancies.  The  energy  of  the  In-X  center  is  between  that  of 
aluminum  and  indium  which  has  suggested  Al-in  distant  pairs. 

During  the  course  of  this  study,  a  number  of  experiments  were  per¬ 
formed  to  pin  down  the  identification  of  the  X-center.  This  work 


7R.  Baron,  M.H.  Young,  J.K.  Neeland,  and  O.J.  Marsh,  Appl .  Phys . 
Lett.  30,  594  (1977). 

®M.W.  Scott,  Appl .  Phys ■  Lett .  32 ,  540  (1978). 

9V.  Swaminathan,  J.E.  Lang,  P.M.  Heminger ,  and  S.R.  Smith,  Appl . 
Phys.  Lett.  35 ,  184  (1979). 

10M.C.  Ohmer  and  J.E.  Lang,  Appl.  Phys.  Lett.  34,  750  (1979). 
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is  described  in  Appendices  C  and  D  which  are  based  on  articles 
submitted  for  publication.  The  results  are  reviewed  here.  Two 

main  conclusions  have  come  out  of  this  work  on  the  X-center.  The 
first  is  that: 

•  There  is  a  class  of  X-center  defects,  with  one  X-center 

associated  with  each  of  the  group  IIIA  elements  in  silicon. 

With  the  indent if ication  of  the  indium-X  center  as  being  an  indium- 
related  complex  with  an  energy  level  about  three-fourths  that  of 
the  subsititut ional  indium,  we  have  been  searching  for  similar 
defects  in  silicon  doped  with  the  other  group  IIIA  acceptors. 

Optical  absorption  has  shown  similar  defects  in  boron-doped  11 

8  11  9 
aluminum-doped,  and  gallium-doped,  silicon.  The  presence  of 

these  defects  has  helped  to  explain  discrepancies  sometimes  ob¬ 
served  in  the  activation  energies  measured  for  these  dopants  by 
Hall  techniques  and  by  optical  absorption.8,11  The  indium-X  level 
can  also  be  seen  in  deep-level  transient  spectroscopy  data.  This 
is  reported  in  Appendix  B.  A  similar  defect  has  now  been  identified 
in  thallium-doped  silicon,  deep-level  transient  spectra  (see 
Appendix  D),  completing  the  set  for  the  group  IIIA  elements. 

These  X-centers  are  all  electrically  active  acceptors  having 
energies  ^80%  of  the  corresponding  group  IIIA  substitutional 
element's  energy.  The  annealing  behavior  for  the  optically  observ¬ 
ed  centers  has  been  measured.  None  of  the  X-centers  anneals  out 
up  to  the  melting  point  of  silicon.  The  concentrations  decrease 
with  higher  annealing  temperature  or  increase  with  lower  anealing 
temperatures.  In  all  cases  this  can  be  done  cyclically.  The  aneal¬ 
ing  kinetics  for  a  number  of  samples  with  different  doping  levels 
can  be  fit  to  that  expected  for  the  reaction 

11W.  Scott  and  C.E.  Jones,  J.  Appl .  Phys .  50,  7258  (1979). 
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(1) 


A  +  C  |  X 


where  A  is  the  group  IIIA  acceptor,  C  is  carbon,  and  K  is  the  equil¬ 
ibrium  constant  equal  to 

K  =  BS  =  K  e-AS/KeAH/KT  (2) 

[A]  [C] 

where  AS  is  the  change  in  the  entropy  and  AH  is  the  change  in  en¬ 
thalpy  for  the  reaction  in  Equation  (1).  These  data  are  shown  in 
Figure  8.  The  AH  term  is  approximately  equal  to  the  X-center 


TEMPERATURE  (°C) 


Figure  8.  Normalized  equilibrium  concentrations  of  X-centers 
as  a  function  of  temperature.  The  quantity  (X)/(A) 
(C)  is  equal  to  the  equilibrium  constant  for  the 
reactions  A  +  C  +  X, 
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binding  energy  and  these  have  been  measured  as  being  -0.03  eV  for 
boron-X,  0.28  eV  for  aluminum-X,  0.40  eV  for  gallium-X,  and  0.70  eV 
for  indium-X.  The  properties  of  the  X-centers  are  listed  in  Table  1. 


TABLE  1.  X-CENTER  PROPERTIES 


Dopant 

E(  X) 
(meV) 

E(ACC) 

Substitutional 

Dopant 

(meV) 

E(X) /E( ACC) 

[Xj/[ACC] 
Pulled  Si 

ea 

(eV) 

B 

37.1 

44.3 

0.84 

io-2 

0.  03 

A1 

56.3 

68.5 

0.82 

10_2-10-3 

0.28 

Ga 

57.0 

72.7 

0.79 

10"4-10-5 

0.40 

In 

112.8 

156.0 

0.72 

10_2-10~4 

0.70 

T1 

180.  0 

246.0 

0.73 

io-2 

- 

It  is  the  similar  X  to  substitutional  acceptor  energies,  the 
similar  annealing  properties,  and  the  defect  binding  energies 
which  go  in  a  smooth  progression  which  have  shown  that  these  de¬ 
fects  are  members  of  a  class  of  complexes  of  the  same  type. 

The  second  result  of  the  X-center  studies  is  that: 

•  The  model  which  best  fits  the  X-center  data  is  that 
of  the  nearest-neighbor  carbon-acceptor  pair. 

The  In-X-center  concentration  has  been  found  to  vary  linearly  with 
both  indium  and  carbon.  No  dependences  can  be  found  on  oxygen  or 
boron-,  and  aluminum  and  thallium  were  not  detected.  These  data 
are  presented  in  Appendix  D.  The  aluminum-X  stress  symmetry  has 
been  determined  to  be  that  of  a  [ill]  trigonal  defect  (see  Appen- 
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dix  C) .  This  rules  out  both  substitutional  single-atom  defects  and 
distant  pairs.  The  high-temperature  stability  and  the  cyclic 
annealability  of  the  X-centers  is  not  consistent  with  interstitial 
or  vacancy-related  defects.  All  of  the  known  interstitial  or 
vacancy-related  defects  anneal  within  a  few  hundred  degrees  of 
room  temperature  and  they  cannot  be  cyclically  annealed.  The  high- 
temperature  stability  and  cyclic  annealing  do  fit  a  model  of  a 
complex  formed  with  substitutional  atoms.  The  [ill]  symmetry  and 
the  dependence  on  a  single  carbon  and  a  single  indium  fit  the 
nearest -neighbor  substitutional  carbon-acceptor  pair  model.  This 
model  is  shown  in  Figure  9. 


Figure  9.  Model  of  the  X-center  as  being  a 
nearest-neighbor  substitutional 
carbon-acceptor  pair.  Sites  1, 
2,3,  and  4  are  different  config¬ 
urations  for  the  carbon. 
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SECTION  4 
TRAPPING  STUDIES 

Deep-level  transient  spectroscopy12’ 13  has  been  used  to  investigate 
the  trapping  parameters  of  defects  in  indium-doped  silicon.  The 
In-X-center  substitutional  indium  and  two  deeper  hole  traps  at 
Ev  +  °-31  eV  and  Ev  +  0-45  eV  were  studied.  The  energy  levels, 
concentrations  emission  rates,  and  capture  cross  sections  were 
determined  for  these  defects.  These  data  were  formalized  and  sub¬ 
mitted  for  publication  and  are  presented  in  Appendix  B.  This  work 
is  briefly  reviewed  below. 

Deep-level  transient  spectroscopy  (DLTS )  measures  the  capacitance 
of  a  reverse-biased  diode.  This  capacitance  depends  on  the  charge 
in  the  depletion  region  of  the  diode.  In  reverse  bias,  most  of  the 
traps  in  this  region  are  empty.  Short  voltages  pulses  are  used  to 
bring  carriers  into  this  region,  where  some  of  the  carriers  are 
trapped.  After  the  pulse,  the  trapped  carriers  are  thermally  emit¬ 
ted,  producing  a  transient  in  the  measured  capacitance.  Electronic 
processing  is  used  to  decompose  the  transient  into  a  spectral  form, 
with  a  peak  for  each  defect  displayed  as  a  function  of  temperature 
and  emission  rate.  The  size  of  the  peak  can  be  related  to  the  trap 
density,  the  change  in  the  emission  rate  as  a  function  of  tempera¬ 
ture  gives  the  trap  activation  energy,  and  the  peak  amplitude  as 
a  function  of  the  filling  pulse  time  is  related  to  the  trap  capture 
cross  section.  A  DLTS  spectrum  for  indium-doped  silicon  is  shown 
in  Figure  10. 


12 

D.V.  Lang,  J.  Appl .  Phys .  45,  3014,  3023  (1974). 

13 

G.L.  Miller,  D.V.  Lang,  and  L.C.  Kimmerling,  in  Ann.  R<^v  of 
Nat.  Sci.  1,  377  (1977).  - - - r 
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Figure  10.  DLTS  spectrum  for  indium-doped  silicon. 


Electric  fields  have  been  found  to  lower  the  thermal  activation 

energies  and  thus  increase  the  emission  rates  for  both  the  In-X 

center  and  substitutional  indium.  These  effects  were  measured 

experimentally  as  a  function  of  electric  field  and  temperature. 

The  data  could  be  fit  by  a  theoretical  model  developed  by 
14 

Hartke  but  the  theoretical  correction  to  zero  field  gives  DLTS 
trap  activation  energies  lower  than  optically  measured  energies. 

The  nitrogen  temperature  capture  coefficients  were  measured  as 

—  7  S  —  7  3 

being  7.6  x  10  cm  /sec  for  indium  and  1.5  x  10  cm  /sec  for 

the  In-X  center.  The  temperature  dependence  for  these  coefficients 

was  estimated  as  being  0.03  ±  0.02  eV  for  indium  and  0.006  ±  0.005 

eV  for  the  In-X  center.  The  capture  coefficients  for  the  hole 


14J.L.  Hartke,  J.  Appl.  Phys .  39,  4871  (1968). 
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traps  at  E  +  0.31  eV  and  E  +  0.45  eV  were  estimated  as  being 
—9  v3  v  _g  3 

2  x  10  cm  /sec  and  1.2  x  10  cm  /sec,  respectively. 


Figure  11  shows  the  emission  rates  as  a  function  of  temperature 

for  these  traps  and  their  effective  thermal  activation  energies. 

The  indium  and  indium-X-center  data  are  also  shown  corrected  to 

14 

zero  field  by  Hartke's  theory. 


I/T  X  1000 


25  12-5  8.33  6  25  50 


Figure  11.  The  emission  rate  of  traps  observed  in  indium-doped 
silicon  as  a  function  of  temperature.  The  indium 
and  indium-X  center  data  are  also  shown  corrected  to 
zero  field  by  the  theory  developed  by  Hartke.14 
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SECTION  5 

ACCOMPLISHMENTS  AND  CONCLUSIONS 


Significant  progress  was  made  in  developing  the  solution-growth 
technique,  in  identifying  and  eliminating  the  In-X  center,  and  in 
measuring  the  trapping  parameters  of  defects  in  indium-doped 
silicon . 

Solution  growths  provided  a  means  of  obtaining  high  indium  con¬ 
centrations  with  good  uniformity  and  with  low-boron  and  X-center 
concentrations.  The  technique  can  be  used  for  volume  production 
adequate  for  projected  detector  needs.  The  crystal  quality  im¬ 
proved  greatly  over  the  period  of  the  contract  but  2-inch-diameter 
wafers  still  show  some  voids  or  polycrystalline  areas. 

The  X-centers  were  identified  as  a  class  of  defects  present  in 
silicon  doped  with  any  of  the  group  IIIA  acceptors.  The  energy 
levels,  optical  spectrum,  annealing  properties,  and  binding 
energies  were  determined  for  most  of  these  centers.  The  X-centers 
were  found  to  depend  on  carbon  and  the  best  model  for  the  defect 
was  identified  as  a  nearest-neighbor  carbon-acceptor  pair.  The 
X-center  can  be  eliminated  by  using  low-carbon  silicon  starting 
materials  in  the  solution  growth  technique. 

The  trapping  parameters  for  substitutional  indium,  the  indium-X 
center,  and  two  deeper  indium-related  traps  were  measured  by 
deep-level  transient  spectroscopy  techniques  near  the  detector 
operating  temperatures.  The  electric  field  effects  were  determined 
for  indium  and  the  indium-X  centers,  and  capture  coefficients 
were  measured  for  all  four  of  the  traps. 
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The  specific  accomplishments  of  the  study  effort  are  listed  as 
follows : 


•  The  solubility  of  indium  in  silicon  was  determined  as 
a  function  of  temperature. 

•  Indium-doped  silicon  was  grown  by  solution  growth  at  the 

18  3 

maximum  of  the  indium  solubility  curve  at  1.6  x  10  In/cm 

8  3 

with  less  than  10  electrically  active  boron  atoms/cm  . 

13  3 

•  The  X-centers  were  eliminated  (less  than  5  x  10  In-X/cm 

g 

by  optical  measurements,  less  than  10  electrically 
active  by  Hall). 

•  X-centers  related  to  boron,  aluminum,  gallium,  indium, 
and  thallium  were  identified. 

•  The  X-center  concentrations  as  a  function  of  acceptor 
concentration,  carbon  concentration,  and  annealing 
temperature  were  determined. 

•  X-center  binding  energies  of  0.03,  0.28,  0.40,  and 
0.70  eV  were  measured  for  the  boron,  aluminum,  gallium, 
and  indium-X  centers,  respectively. 

e  The  nearest-neighbor  carbon-acceptor  pair  model  was 
identified  as  the  best  model  for  the  X-center  complex. 

•  The  trapping  parameters  for  indium,  the  indium-X  center, 

and  two  indium-related  hole  traps  at  Ev  +  0.31  eV  and 

E  +  0.45  eV  were  determined, 
v 
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APPENDIX  A 

SOLUTION  GROWTH  OF  INDIUM-DOPED  SILICON* 


ABSTRACT 

Indium-doped  silicon  has  been  grown  from  indium-rich  solutions 

using  a  gradient-transport  solution  growth  process.  The  growth 

temperatures  were  varied  from  950°  to  1300°C  to  determine  the 

solubility  limits  of  indium  in  silicon.  The  maximum  indium  con- 

cent ration  obtained  was  1.6  x  10  /cm  at  a  growth  temperature  of 

1300°C  but  indications  are  that  the  maximum  solubility  is  2.5  x 

10  /cm  .  The  oxygen  and  carbon  concentrations  increased  with  the 

increased  growth  temperatures,  suggesting  solubility-limited  values 

The  shallower,  indium-X  defect  also  increased  with  growth  tempera- 

13  3 

ture,  but  could  be  reduced  below  2  x  10  /cm  by  removing  carbon 
from  the  crystals.  The  uniformity  in  indium  doping  was  typically 
±5%  across  the  crystal  diameter. 

INTRODUCTION 

Indium-doped  silicon  is  being  investigated  as  an  extrinsic  photo¬ 
conductor  material  for  use  in  the  3-  to  5-ym  infrared  region.  A 
number  of  problems  have  arisen  with  this  material,  notably  with 


*This  work  was  sponsored  in  part  by  the  Defense  Advanced  Research 
Projects  Agency  under  Order  No.  3211,  monitored  by  NV&EOL  under 
Contract  No.  DAAK-77-C-0194 .  The  authors  are  Walter  M.  Scott, 

Colin  E.  Jones,  and  Robert  J.  Hager,  Honeywell  Corporate  Material 
Sciences  Center,  Bloomington,  MN.  The  paper  was  presented  at  the 
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the  low  indium  concentration  obtained  by  standard  growth  methods 

and  the  presence  of  a  shallower  defect  level  associated  with  the 

indium.  Czochralski-grown  crystals  have  been  limited  to  a  maximum 
•  •  17  3 

indium  concentration  of  3.8  x  10  /cm  before  microprecipitates 

A-l 

are  observed  in  the  crystals.  The  concentration  of  indium  in 
float-zone-grown  crystals  is  generally  even  lower  than  this  be¬ 
cause  of  the  low  segregation  coefficient  of  indium. A-2  The  limits 
to  the  indium  solubility  in  silicon  have  not  been  previously  re¬ 
ported  in  the  literature  to  determine  what  improvement  could  be 
made  in  the  doping  levels  in  silicon  crystals. 

A  shallower  defect  center  associated  with  both  indium  and  car- 
A-3  A-4 

bon,  ’  and  labeled  as  indium-X,  has  been  observed  in  indium- 

A-5 

doped  silicon.  Baron  et  al .  observed  this  defect  in  Hall  co¬ 
efficient  measurements  at  an  energy  of  about  0.11  eV.  Subsequently, 
A-6 

Scott  observed  the  infrared  spectrum  of  an  effective-mass-like 
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E.L.  Kern,  R.  Baron,  R.H.  Walker,  D.J.  O'Connor,  and  O.J. 

Marsh,  J.  Electron.  Mater.  4,  1249  (1975). 

A- 2 

W.R.  Runyan,  Silicon  Semiconductor  Technology,  (McGraw  Hill, 

New  York,  1965 ) . 

A- 3 

C.E.  Jones,  D.  Schafer,  W.  Scott,  and  R.J.  Hager,  Bull.  Am. 
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defect  associated  with  indium  at  an  energy  of  0.1128  eV.  The  origin 
of  the  defect  is  uncertain  at  this  time,  but  it  is  known  to  exist 
in  both  Czochralski  and  float-zone-grown  crystals . 

Solution  growth  techniques  such  as  gradient-transport  solution 

growth  and  liquid-phase  epitaxy  (LPE)  have  been  used  to  grow 

group  III-V  compounds  for  many  years,  but  these  techniques  have 

A-7 

not  been  extensively  applied  to  silicon.  Pfann  described  the 

use  of  temperature-gradient  zone  melting  as  a  means  of  producing 

large-area  p-n  junctions  in  silicon,  and,  more  recently,  Cline 
A-8 

and  Anthony  have  explored  the  migration  of  liquid-metal  drop¬ 

lets  in  silicon.  An  important  feature  of  these  solution  growth 
procedures  is  that  the  crystals  are  grown  below  the  melting  point 
of  silicon,  so  the  retrograde  solubility  of  the  impurities  can  be 
used  to  obtain  very  heavily  doped  crystals. 

In  this  paper  we  report  the  results  of  gradient-transport  solution- 

growth  of  silicon  from  indium  solutions.  The  growths  were  done  as 

a  function  of  temperature  to  determine  the  solubility  limits  of 

indium  in  silicon.  Growth  rates  of  about  4  mm/day  have  been  achieved 

18  3 

and  crystals  doped  up  to  1.6  x  10  /cm  with  indium  have  been  pro¬ 
duced  by  this  technique.  Single  crystals  in  excess  of  2.5  cm  diam¬ 
eter  have  been  produced  with  a  doping  uniformity  of  ±5%.  Indium-X 

13  3 

concentrations  <5  x  10  /cm  have  been  obtained  in  crystals  with 

„  iri17  T  .  3 

7  x  10  In/ cm  . 


A  7W.G.  Pfann,  Solid  State  Physics,  Volume  4,  F.  Seitz  and  D. 
Turnbull,  Eds.  (Academic  Press,  N.Y. ,  1957)  pp .  424. 

A-8H.E.  Cline  and  T.R.  Anthony,  J.  Appl .  Phys .  48,  2196  (1977), 
and  the  references  therein. 
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CRYSTAL  GROWTH  PROCEDURES 


Solubility  Limits 

A  schematic  illustration  of  the  growth  apparatus  used  to  determine 
the  solubility  limit  of  indium  is  shown  in  Figure  A-l.  It  consists 
of  a  flat-bottomed  quartz  tube  in  which  was  placed  a  single-crystal 
silicon  seed,  high-purity  (6N)  indium,  and  a  high-purity  silicon 
source  in  the  relative  positions  shown  in  Figure  A-l.  A  tubular 
quartz  spacer  was  used  to  hold  the  silicon  seed  in  place  during 
the  growth  and  prevent  it  from  buoying  to  the  surface  of  the  liquid 
indium.  The  apparatus  was  placed  in  a  furnace  with  uniform  temper¬ 
ature  profile  and  a  temperature  difference  established  between  the 
silicon  source  and  the  substrate  by  a  flow  of  air  directed  at  the 
flat  bottom  of  the  ampoule.  The  temperature  of  the  silicon  seed  is 
maintained  approximately  50°C  colder  than  the  source  temperature 
throughout  the  growth  process . 
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Figure  A-l.  Schematic  of  the  solution-growtn  apparatus. 
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The  temperature  difference  between  source  and  substrate  was 
determined  by  calibrating  the  system  using  a  Pt-Pt  (13%  Rh)  thermo¬ 
couple  .  The  capsule  was  loaded  with  the  usual  amount  of  indium 
and  a  standard  substrate.  A  quartz  tube  containing  the  thermo¬ 
couple  was  inserted  into  the  indium  to  be  in  contact  with  the  sub¬ 
strate.  The  apparatus  was  brought  up  to  temperature,  and  the 
temperature  profile  in  the  indium  measured.  The  profile  was  mea¬ 
sured  with  the  gas  flow  rate  a  parameter  in  order  to  be  able  to 
set  any  desired  temperature  difference  across  the  indium  column. 

The  repeatability  of  the  temperature  setting  was  estaimated  at 
about  ±5°C  from  these  calibration  measurements. 

Large-Area  Growth 

A- 9 

The  growth  configuration  was  modified  to  grow  larger-area 
crystals  (up  to  2-inch  diameter).  The  length  of  indium  column 
was  reduced  to  about  0.05  cm  by  containing  the  solution  in  a  well 
formed  in  the  silicon  substrate.  Other  details  of  the  growth  re¬ 
mained  basically  the  same  as  shown  in  Figure  A-l. 

RESULTS  AND  DISCUSSION 

Electrical  Evaluation 


All  the  crystals  were  characterized  by  Hall  coefficient  versus 
temperature  measurements  on  at  least  one  sample  from  each  growth 
taken  from  the  single-crystal  portion  of  the  growth.  A  major 
problem  in  the  electrical  transport  measurements  of  deep  im¬ 
purities  in  silicon  is  in  determining  the  total  concentration  of 

A-^c.e.  Jones,  D.E.  Schafer,  M.W.  Scott,  and  R.J.  Hager,  "Studies 
of  Indium-Doped  Silicon,"  Semiannual  Report,  NVL  Contract  No. 
DAAK70-77-C-0194 ,  31  July  1978. 
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the  impurities.  Hall  measurements  provide  a  measure  of  the  free- 
carrier  concentration,  but  if  all  the  impurities  are  not  ionized 
at  the  temperature  of  measurement,  then  a  specific  model  must  be 
adopted  before  the  total  impurity  concentration  can  be  determined 

The  free-carrier  concentration,  p,  was  determined  from  the  Hall 
coefficient,  Rjj,  from  the  expression 


where  the  scattering  factor,  r,  was  assumed  to  be  unity.  The 
carrier  concentration  is  also  determined  from  the  charge  neutral¬ 
ity  condition 

P  +  nd  =  NA 


where  Np  is  the  concentration  of  ionized  donors,  and  N~  is  the 
concentration  of  ionized  acceptors.  By  assuming  a  low  concentra¬ 
tion  of  compensating  donors  and  a  high  concentration  of  indium 
relative  to  any  other  acceptor,  this  can  be  rewritten  as 


N 


P  = 


In 


1  +  g  exp  — 


E.  -  E, 


KT 


(A-l) 


where  N^n  is  the  total  indium  concentration,  g  is  the  ground- 

state  degeneracy,  E.  is  the  indium  ionization  energy,  and  E„  is 

i  F 

the  Fermi  energy.  In  this  expression  the  effects  of  the  indium- 
excited  states  have  been  ignored,  and  only  the  ground  state 
considered . 
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The  hole  concentration  is  also  determined  from  the  usual  expres¬ 
sion  for  nondegenerate  bands 


N 


P  = 


v 


EF 

1  exp^ 


(A-2  ) 


where  N  is  the  effective  density  of  states  in  the  valence  band, 
v 

Combining  Equations  (A-l)  and  (A-2)  gives  the  hole  concentration 
as  a  function  of  the  total  indium  concentration. 

The  values  of  the  various  parameters  used  in  the  calculation  are 
the  following: 

N  =  4.8296  x  1015  (m*  /m  T)3/2 
v  h'  o 

m*  =  0.81  at  300K  (ref.  A- 10) 
m 

E.  =  0.155  eV 

l 

g  =4 

The  properties  of  the  various  crystals  are  summarized  in  Table 
A-l.  The  300K  electrical  properties  are  listed  along  with  the 
growth  temperature  and  the  calculated  indium  concentration.  The 
boron  concentration  observed  in  the  crystals  is  also  listed  in 
Table  A-l,  and  was  obtained  from  the  exhaustion  plateau  observed 
in  the  Hall  coefficient  versus  reciprocal  temperature  plot .  The 
cases  where  no  boron  concentration  is  listed  are  those  where  the 
Hall  coefficient  versus  reciprocal  temperature  curve  did  not  have 
a  boron  exhaustion  plateau. 

A-10H.D.  Barber,  Solid  State  Elect.  10,  1039  (1967). 
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TABLE  A-l.  SUMMARY  OF  ELECTRICAL  EVALUATIONS 


Crystal 

No. 

Growth 

Temp 

<°C) 

Rtt(  300  ) 

H  3 

(cnr /C) 

p ( 300 ) 
(fi-cm) 

Indium 

/  -3n 

(cm  ) 

Boron 

(cm-3) 

1 

946 

620.0 

2.17 

1.8  x  1016 

2 

1005 

281.0 

1.28 

o 

o 

X 

I—1 

o 

I—1 

<S) 

3 

1008 

317.0 

1.03 

5.2  X  1016 

14 

1.0  x  10x* 

4 

1056 

169.0 

0.60 

17 

1.5  x  10 

2.0  x  1015 

5 

1056 

198.0 

0.74 

17 

1.1  x  10  ' 

6 

1056 

213.0 

0.75 

1.0  x  1017 

14 

5.0  x  10x* 

7 

1100 

122.0 

0.44 

17 

2.8  x  10  ' 

1.5  x  1015 

8 

1150 

87.5 

0.36 

17 

5.1  x  10  ' 

1.0  x  1015 

9 

1150 

79.0 

0.32 

17 

5.2  x  10x/ 

10 

1150 

75.0 

0.32 

6.8  x  1017 

11 

1200 

62.4 

0.24 

1.0  x  1018 

6.3  x  1015 

12 

1200 

57.0 

0.26 

1.1  x  1018 

13 

1250 

52.5 

0.26 

1.4  x  1018 

14 

1300 

45.5 

0.24 

1.6  x  1018 

15 

1300 

45.0 

0.21 

1.7  x  1018 
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A-ll 

Recently,  Sclar  published  an  analysis  relating  the  300K 

resistivity  of  silicon  to  the  concentration  of  deep  impurities 
in  the  crystal.  The  calculation  takes  into  account  the  variation 
in  mobility  with  free-carrier  concentration  and  the  usual  mass- 
action  equation  to  calculate  the  free-carrier  concentration.  His 
curve  of  300K  resistivity  versus  indium  concentration  is  repro¬ 
duced  in  Figure  A-2  along  with  the  values  obtained  here.  As  is 
obvious  from  the  figure,  the  procedure  used  here  agrees  well  with 
Sclar 's  calculated  curve. 


Figure  A-2.  Room  temperature  resistivity  of  indium-doped 
silicon  as  a  function  of  the  indium  concen¬ 
tration.  The  circles  are  the  measured  values 
and  the  solid  curve  is  taken  from  the  calcu¬ 
lations  of  Sclar  (ref.  A-ll). 


A  1:LN.  Sclar,  IEEE  Trans.  Electron  Devices  ED-24,  709  (1977). 
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Solubility  Limit 


The  indium  concentration  in  the  solution-grown  crystals  is  shown 

in  Figure  A-3  as  a  function  of  the  growth  temperature.  The  maxi- 

1  ft 

mum  value  of  indium  obtained  in  our  growths  has  been  1.6  x  10  / 

cm2  at  1300°C.  Higher  growth  temperatures  were  not  attempted 
because  of  furnace  limitations  and  problems  with  the  quartz 
ampoules  at  higher  temperatures . 


Figure  A-3.  Measured  indium  concentrations 
in  solution-grown  crystals  at 
various  growth  temperatures. 


A-12 

Hall  has  suggested  that  the  distribution  coefficient,  k,  of 

a  solute  having  a  low  solid  solubility  should  have  a  temperature 
dependence  given  by  an  equation  of  the  type  log  k  =  B  -  A/T,  where 


A  12R.N.  Hall,  J.  Phys.  Chem.  Solids  3,  63  (1957) 
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A  and  B  are  constants,  and  T  is  the  absolute  temperature.  The 
values  for  distribution  coefficient  which  we  have  obtained  are 
shown  in  Figure  A-4,  and  as  is  evident,  they  correspond  to  this 

type  of  relationship.  The  distribution  coefficient  for  indium  at 

—4 

the  silicon  melting  point,  k  ,  is  estimated  to  be  2.5  x  10  from 

°  -4 

our  data.  This  is  significantly  lower  than  the  value  of  4  x  10 

usually  used  for  indium.  It  was  pointed  out  that  this  simple 


Figure  A-4.  Distribution  coefficient,  k,  of  indium  in 
silicon  as  a  function  of  reciprocal  temp¬ 
erature.  The  distribution  coefficient  at 
the  melting  point,  kQ,  is  2.5  x  10  4 . 


A-13F.A.  Trumbore,  Bell.  Syst .  Tech.  J.  39,  205  (1960). 
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equation  does  not  necessarily  describe  the  temperature  dependence 
of  the  distribution  coefficient  but  empirically  it  has  been  found 
for  a  number  of  impurities  in  silicon  and  germanium.  The  fact  that 
most  of  the  data  in  Figure  A-4  lie  on  a  straight  line  would  suggest 
that  this  relation  applies  to  indium  in  silicon.  In  this  case  the 
crystals  at  the  upper  two  temperatures  may  have  been  undersatur¬ 
ated.  The  solubility  maximum  estimated  from  the  straight  line  is 

183  in  o 

about  2.5  x  10  /cm  ,  rather  than  the  1.6  x  10  /cm  which  we  have 
obtained . 

Doping  Uniformity 

A  major  concern  in  producing  material  for  focal  planes  is  the 
uniformity  of  doping.  This  has  been  determined  on  solution-grown 
crystals  using  spreading  resistance  measurements,  and  the  results 
are  shown  in  Figure  A-5  along  with  the  results  from  Czochralski 
and  float-zone  indium-doped  material.  The  solution-grown  crystal 
is  uniform  to  ±5%,  or  better.  The  Czochralski  growth  shows  a 


iioo 

S 

<2 

* 

if 


Figure  A-5.  Spreading  resistance  across  a  solution- 
grown  crystal  compared  with  spreading 
resistance  across  Czochralski  and  float- 
zone-grown  crystals. 
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drift  in  doping  level  from  the  center  to  the  edge  along  with  some 
striations.  Total  variations  across  the  entire  crystal  are  about 
±20%.  Larger  striations  are  also  seen  in  the  float— zone  material, 
with  a  total  variation  from  center  to  edge  well  in  excess  of  ±20%. 


Growth  Rates 


The  growth  rate  can  be  computed  from  the  standard  expression: 


J 


3C 


-D 


Si 


9x 


substrate 


2 

where  J  is  the  particle  current  in  moles/cm  sec,  and  C  ^  is  the 
silicon  concentration  in  moles/cm^.  The  growth  rate  at  the  sub¬ 
strate  is  given  by 


V(cm/sec)  =  J - 

pSi 

where  is  the  molar  weight  of  silicon,  and  is  the  density 

of  the  silicon  crystal.  In  steady  state  the  growth  rate  can  be 
written 


...  .  ,  “si  <Co  -  Cl> 

V(cm/sec)  =  D -  - : - 

pSi  X 


since  the  concentration  gradient  will  be  constant  throughout  the 
solution.  To  characterize  the  growth  process  completely,  we  need 
to  know  only  two  parameters:  (a)  the  diffusion  coefficient,  D, 
and  (b)  the  silicon-indium  liquidus  (to  obtain  Cg^).  The  diffu¬ 
sion  coefficient  in  liquid  indium  has  been  measured  up  to  a 
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u 


temperature  of  1000°C,  and  is  given  by  the  expression^  ^ 
D  =  3  x  10-4exp  -  hS.  kcal/mole  om2/sec 


where  R  is  the  gas  constant.  The  silicon-indium  liquidus  values 
used  were  obtained  from  the  data  of  Thurmond  and  Kowalchick. A-15 

The  predicted  growth  rate  in  units  of  mm/day  is  shown  in  Figure 
A-6  as  a  function  of  the  growth  temperature.  For  a  given  temper- 


Figure  A-6.  Growth  rates  for  diffusion  through  a  3-cm  column 
of  indium  with  a  50°C  temperature  difference 
from  source  to  substrate.  The  solid  line  is  cal¬ 
culated  from  a  one-dimensional  diffusion  model. 


The  data  for  diffusion  in  liquids  are  summarized  in:  N.H. 
Nachtreib,  Ady .  Phvs .  16 ,  309  (1967). 

A-15 

C.D.  Thurmond  and  M.  Kowalchick.  Bell  Svst .  Tech.  J.  39 
169  (1960)  ’ 
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ature  gradient,  the  growth  rate  increases  dramatically  with  in¬ 
creasing  growth  temperature  up  to  about  1300°C  because  of  the 
increased  solubility  of  silicon  in  indium  at  the  higher  temper¬ 
ature.  The  slope  of  the  liquidus,  has  a  maximum  near  1300°C, 

accounting  for  the  maximum  in  growth  rate.  The  maximum  rate  is 
expected  to  be  in  the  10  to  20-mm/day  range  under  these  growth 
conditions . 


The  measured  growth  rates  are  also  shown  in  Figure  A-6 .  In  general, 
the  measured  values  agree  very  well  with  calculated  values,  in¬ 
dicating  that  reasonable  values  of  the  parameters  were  used  in  the 
calculation.  The  errors  assigned  to  the  growth  rates  result  from 
an  uncertainty  in  the  degree  of  melt-back  of  the  substrate  and 
also  from  an  uncertainty  in  the  amount  of  growth.  On  cooldown, 
material  would  precipitate  on  the  solution-grown  crystal,  giving 
an  uncertainty  as  to  the  length  of  the  crystal .  Some  error  also 
occurs  in  determining  the  length  of  time  of  the  growth,  since  it 
requires  some  time  for  steady  state  to  be  reached.  The  growths 
were  typically  about  10  days,  so  this  is  judged  to  be  only  a 
small  error.  This  uncertainty  in  growth  rate  is  typically  ±20% 
so  still  provides  a  sufficiently  accurate  comparison  with  the 
simple  one-dimensional  diffusion  model. 

Optical  Evaluation 

Optical  absorption  measurements  were  used  to  monitor  the  concen¬ 
tration  of  oxygen,  carbon,  boron,  indium  and  the  shallower  indium- 
X  defect  in  the  crystals,  and  also  determine  if  other,  shallower 
impurities  are  introduced  into  the  crystals  by  the  solution 
growth  process.  The  absorption  spectrum  of  one  of  our  solution- 
grown  indium-doped  silicon  crystals  in  the  400-  to  2000-cm  ^ 
region  is  shown  in  Figure  A-7 .  The  principal  indium  lines  are 
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marked  as  lines  1,  2,  3,  and  4,  to  be  consistent  with  the  number- 

A  “| 

ing  system  used  by  Onton  et  al. "  for  acceptors  in  silicon. 

The  largest  absorption  line  observed  at  1134  cm-1  is  due  to 
oxygen.  The  broad  absorption  bands  in  the  400-  to  1000-cm-1  re¬ 
gion  are  lattice  absorptions,  and  the  broad  absorption  starting 

at  1200  cm  is  due  to  ionization  of  the  indium  centers.  The 
region  of  interest  for  monitoring  the  indium-X  defect  is  between 
700  and  1000  cm"  ,  where  there  is  appreciable  lattice  absorption. 
The  lattice  absorption  bands  were  subtracted  from  the  absorption 
spectrum  as  described  earlier  to  observe  these  defect  absorption 
lines.  ■L/ 


Figure  A-7.  Infrared  absorption  spectrum  of  silicon  containing 
1  x  1017  indium/cm3.  The  major  indium  lines  are 
numbered  1  through  4,  and  the  absorption  lines  due 
to  oxygen  and  carbon  are  also  indicated. 


A- 16 

A.  Onton,  P.  Fisher, 
(1967). 


and  A.K.  Ramdas,  Phys .  Rev.  163,  686 


A-17 

M.W.  Scott,  C.E.  Jones,  and  R.J.  Hager,  "Studies  of  Indium- 
Doped  Silicon,"  Semiannual  Report,  NVL  Contract  No  DAAK70- 
77-C-0194,  31  Jan  1978. 
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The  shallower  indium-X  defect  was  previously  shown  to  have  an 

effective-mass-like  absorption  spectrum  with  the  principal  ab- 

-1  A -6 

sorption  lines  at  800,  830,  and  873  cm  .  The  integrated  in¬ 
tensity  of  the  strongest  line  at  830  cm  ^  was  used  to  determine 
the  concentration  of  this  defect.  The  concentration  of  the 

indium-X  defect  in  the  as-grown  crystals  is  listed  in  Table  A-2 . 

3  4 

Since  it  is  known  ’  that  the  concentration  of  the  indium-X 
defects  varies  with  annealing  temperature,  the  values  in  Table 
A-2  only  represent  the  as-grown  conditions.  An  important  result 
to  note  in  Table  A-2  is  the  fact  that  the  indium-X  concentration 
can  be  made  negligibly  small  in  crystals  where  the  carbon  con¬ 
centration  is  kept  low,  as  in  crystal  9. 


TABLE  A-2.  SUMMARY  OF  CRYSTAL  PROPERTIES  FROM  OPTICAL  DATA 


Crystal 

Growth 

Temp. 

(°C) 

[in] 

[X] 

[o] 

[c] 

W 

3 

1008 

1.3  x  1016 

<1.0  x  1013 

4.5  x  10^ 

1.3  x  1016 

1.4  x  1014 

4 

1056 

1.7  x  1017 

<4.5  x  1013 

j  8.0  x  1016 

3.7fe  x  1016 

2.35  x  1015 

6 

1056 

1.76  x  1016 

<1.0  x  1013 

7.7  x  1016 

2.35  x  1016 

14 

5.0  x  10 

7 

1100 

17 

4.2  x  10 

14 

4.5  x  10 

17 

1.7  x  10±/ 

17 

1.0  x  10-1-' 

5.5  x  1015 

8 

1150 

17 

6.9  x  10A/ 

14 

3.4  x  10 

17 

2.8  x  10 

17 

1.2  x  10  ' 

1.2  x  1016 

9 

1150 

17 

5.24  x  10 

<9.0  x  1012 

2.13  x  1017 

<3.3  x  1015 

<2.5  x  1013 

13 

1250 

17 

8.0  x  10±/ 

1.9  x  1015 

17 

4.5  x  10 

17 

2.2  x  10  ' 

14 

2.7  x  10 

14 

1300 

17 

9.5  x  lcr' 

2.2  x  1015 

17 

6.2  x  10 

17 

2.0  x  10±/ 

<1.0  x  1013 

16 

1150 

5.5  x  1017 

5.0  x  1013 

5.0  x  1015 

1.0  x  1017 

<5.0  x  1013 
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The  concentrations  of  both  oxygen  and  carbon  were  also  measured 
in  the  crystals  as  a  measure  of  crystal  purity.  The  oxygen  con¬ 
centration  was  determined  from  the  peak  intensity  of  the  1134-cm-^ 
line  at  8K,  and  the  carbon  concentration  was  determined  from  the 
intensity  of  the  606-cm-^  line.  The  conversion  from  peak  absorp¬ 
tion  coefficient  to  doping  level  is  taken  from  the  literature 

A-18 

using  the  room  temperature  values  reported  by  Gaff  et  al .  for 

A- 19 

oxygen  and  by  Newman  and  Willis  for  carbon.  The  absorption 

coefficient  for  these  lines  was  measured  as  a  function  of  temper¬ 
ature  to  determine  the  appropriate  conversion  at  8K.  The  values 

used  for  the  8K  measurements  were  3.1  x  10x  /cm  /cm"  for  oxygen, 
16  3  1 

and  6 . 7  x  10  /cm  /cm  for  carbon. 

The  oxygen  and  carbon  concentrations  are  shown  in  Figure  A-8  as 
a  function  of  the  growth  temperature.  The  oxygen  concentrations 
of  four  Czochralski-grown  crystals  are  also  shown  in  Figure  A-8 
for  comparison.  Both  the  oxygen  and  carbon  concentrations  increase 
with  increasing  growth  temperatures,  suggesting  solubility- 
limited  values.  The  oxygen  concentrations  obtained  here  agree 
reasonably  well  with  the  solubilities  measured  by  Hrostowski  and 
Kaiser. in  the  case  of  carbon,  there  is  significant  scatter 
in  the  data,  which  could  mean  that  all  the  crystals  are  not  at 
the  solubility  limit  at  the  particular  growth  temperature. 


A  18K.  Graff,  E.  Grallath,  S.  Ades,  G.  Goldback,  and  G.  Tolg, 
Solid  State  Elect.  16,  887  (1973). 

A-19R.C.  Newman  and  J.B.  Willis,  J.  Phys.  Chem.  Solids  26,  373 
( 1965  )  . 

A'^9H.J.  Hrostowski  and  R.H.  Kaiser,  Phys.  Rev.  107 ,  969  (1957) 


44 


■s 


IS 

10 


17 

10  ■ 


IS 

io  L 


OXYGEN 


CARBON 


1400  1300  1200  1100  1000 

GROWTH  TEMPERATURE  (*C) 


Figure  A-8.  The  oxygen  and  carbon  concentrations  found 

in  indium- doped  silicon  crystals.  The  oxygen 
concentration  in  solution-grown  crystals  is 
shown  as  X's  and  in  Czochralski  crystals  as 
A's.  The  carbon  concentration  in  solution- 
grown  crystals  is  shown  as  0's. 


SUMMARY 


We  have  demonstrated  that  indium-doped  silicon,  doped  at  the 
solubility  limits  of  indium,  can  be  grown  from  an  indium  solu¬ 
tion  using  the  gradient-transport  process.  Growth  rates  of  4  mm/ 
day  have  been  obtained,  making  the  technique  feasible  as  a  method 
of  producing  heavily  doped  silicon  wafers  in  low-volume  applica¬ 
tions.  The  solubility  limit  of  indium  has,  for  the  first  time, 
been  measured  over  the  950°  to  1300°C  temperature  range.  The  max¬ 
imum  solubility  obtained  in  the  solution-grown  crystals  was 
1.6  x  1018cm3,  but  the  indications  are  that  the  maximum  solubility 
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occurs  at  about  1300°C  with  a  value  of  about  2.5  x  1018/cm3. 
Optical  absorption  measurements  have  verified  that  the  concentra¬ 
tion  of  the  shallower  indium-X  defect  is  lower  in  the  solution- 
grown  crystals  when  the  carbon  concentration  is  kept  low.  Uniform¬ 
ity  of  the  indium  concentration  has  been  shown  to  be  ±5%  across 
the  crystal  diameter,  and  the  boron  concentration  has  been  at 
least  equivalent  to  Czochralski  crystals.  The  solution-growth, 
Czochralski-growth,  and  float-zone  growth  methods  are  compared 
in  Table  A-3. 


TABLE  A-3.  COMPARISON  OF  GROWTH  METHODS 


Crystal 

Solution 

Growth 

Czochralski 

Float 

Zone 

In(cm-3) 

1018 

2  to  3  x  1017 

1  x  1017 

0 ( cm  3 ) 

(Solubility  limit) 

^1018 

<1016 

C(cm-3 ) 

%1016 

>1017 

^1016 

B(cm  3) 

<5  x  1013 

„io14 

„1012 

Ap/p 

±5% 

±20% 

>±20% 
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APPENDIX  B 


DLTS  STUDIES  OF  TRAPPING  PARAMETERS  FOR 
CENTERS  IN  INDIUM- DOPED  SILICON* 


ABSTRACT 


Deep-level  transient  spectroscopy  (DLTS)  has  been  used  to  measure 

the  low-temperature  trapping  parameters  of  defects  in  indium- 

doped  silicon.  Substitutional  indium  at  Ev  +  0.15  eV,  the  indium-X 

center  at  Ev  +  0.11  eV,  and  two  deeper  indium-related  centers  at 

E  +  0.31  eV  and  E  +  0.45  eV  were  studied.  Electric  fields  were 
v  v 

found  to  lower  the  activation  energies  and  increase  the  emission 

rates  for  the  substitutional  indium  and  the  indium-X  center. 

Theoretical  models  including  the  field  effect  on  the  barrier  and 

thermally  assisted  tunneling  were  used  to  fit  the  data.  The  capture 

coefficients  near  liquid  nitrogen  temperature  were  estimated  as 

-9  3 

being  for  substitutional  indium,  C(In)  =  7.6  x  10  cm  /sec  exp 
(+0.031  eV/kT);  for  the  indium-X  center,  C(InX)  =  7.7  x  15  ^  cm^/ 
sec  exp  (+0.006  eV/RT);  for  the  Ev  +  0.31  eV  center,  C(H1)  =  2  x 
10"9  cm^/sec;  and  for  the  Ev  +  0.45  eV  center  C(H2)  =  1.2  x  10  ^ 
cm^/sec . 


*Work  done  under  DARPA  Contract  DAAK70-77-C-0194 .  The  authors  of 
this  paper  are  Colin  E.  Jones  and  Gregory  E.  Johnson,  Honeywell 
Corporate  Material  Sciences  Center,  Bloomington,  MN.  The  paper 
was  submitted  to  J.  Appl .  Phys .  for  publication. 
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INTRODUCTION 


Silicon  doped  with  indium  is  being  studied  as  an  infrared  photo- 
detector  material  in  the  3-5  micron  region.  One  of  the  problems 
with  indium  doping  has  been  the  formation  of  a  shallow  indium- 
related  complex  called  the  X-center.B-1  “  B-4  Holes  thermally 
generated  from  these  centers  increase  the  background  noise  in 
indium-doped  silicon  detectors  unless  the  system  is  cooled  to 
around  45K  to  freeze  out  the  X-centers. 

In  order  to  calculate  the  photodetector  properties  of  indium- 
doped  silicon,  the  trapping  parameters  for  substitutional  indium, 
the  In-X  center,  and  other  defects  common  in  this  material  must 
be  known  for  temperatures  in  the  50K  to  100K  range.  This  paper 
describes  the  use  of  deep-level  transient  spectroscopy  ( DLTS )  to 
measure  the  low-temperature  trapping  parameters  of  indium,  the 
indium-X  center  at  Ev  +  0.11  eV,  and  two  deeper  hole  traps  at 
Ev  +  0.31  eV  and  Ey  +  0.45  eV. 


B  1R.  Baron,  M.H.  Young, 
Lett.  30,  594  (1977). 


J.K.  Neeland,  and  O.J.  Marsh , 


Appl .  Phys . 


B_2 

W.  Scott,  Appl.  Phys.  Lett.  32,  540  (1978). 

B—  3 

R.  Baron,  J.P.  Bankus,  S.D.  Allen,  T.C.  McGill,  M.H.  Young, 

H.  Kimura,  H.V.  Winston,  and  O.J.  Marsh,  J.  Appl.  Phys.  Lett. 
34,  257  (1979)  - — - L - 

B-4 

C.E.  Jones,  D.E .  Schafer,  W.  Scott,  and  R.  Hager,  Bull.  Am. 
Phys.  Soc.  24,  276  (1979)  BL14,  and  IRIS  Speca. 
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EXPERIMENTAL 


B— 5  —  B— 7 

Deep-level  transient  spectroscopy  is  the  technique  used 

to  determine  the  concentrations,  energies,  and  trapping  parameters 
for  defects  in  the  indium-doped  silicon.  The  test  samples  are  n+  p 
diodes  made  by  thermally  diffusing  phosphorus  into  silicon  doped 

1  CZ. 

with  indium  and  boron.  The  indium  concentrations  ranged  from  1(j 
3  13  3 

In/cm  to  10  In/cm  .  For  each  sample,  the  boron  concentration 
was  about  10  times  that  of  the  indium.  The  diodes  are  held  in  re¬ 
verse  bias  which  establishes  a  depletion  region.  Free  carriers 
are  pushed  out  of  this  region  by  the  applied  field  and  the  traps 
are  emptied.  Short  voltage  pulses  are  used  to  collapse  the  deple¬ 
tion  region,  bringing  in  carriers,  some  of  which  are  trapped  at 
defects.  After  the  pulse  goes  away,  these  trapped  carriers  are 
slowly  emitted  thermally.  This  produces  a  transient  in  the  diode 
capacitance  which  is  measured  as  a  function  of  temperature.  Elec¬ 
tronic  processing  is  used  to  separate  the  transient  into  a  compo¬ 
nent  for  each  trap  and  to  display  the  data  in  a  spectral  form.  The 

B— 5  B— 7 

resulting  peak  heights  are  related  to  trap  concentrations.  ~  ’ 

The  peak  positions  give  the  trap  emission  rate  as  a  function  of 

temperature.  The  change  in  emission  rate  as  a  function  of  temper- 

B— 5  B— 7 

ature  is  used  to  determine  trap  activation  energies.  ’  The 

saturation  of  the  peak  heights  as  a  function  of  the  filling  pulse 

B— 6 

times  can  be  used  to  calculate  the  capture  cross  sections. 


B  5D.V.  Lang,  J.  Appl .  Phys .  45,  3014  (1974). 

B_6D.V.  Lang,  J.  Appl.  Phys.  45,  3023  (1974). 

B  -  7 

G.L.  Miller,  D.V.  Lang,  and  L.C.  Kimmerling,  Ann.  Rev, 
of  Mater.  Sci .  7,  (1977). 
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The  DLTS-measured  trap  energies  differ  from  the  trap  energies 
measured  by  Hall  or  optical  techniques  due  to  several  reasons. 

Two  of  the  most  important  are  the  temperature  dependence  of  the 
capture  cross  section  and  the  increase  of  the  trap  emission  rates 
in  an  electric  field.  The  mathematics  for  both  of  these  corrections 
are  reviewed  in  the  next  section.  The  experimental  technique  used 
to  determine  the  field  effect  is  described  here. 

The  electric  field  in  the  diode  depletion  region  goes  from  a  maxi¬ 
mum  at  the  junction  to  zero  at  the  end  of  the  depletion  region. 

The  fields  can  be  calculated  as  a  function  of  position  from 
Poisson's  equation 

l|ix)  „  q|x)  (B-l) 

where  F(x)  is  the  electric  field  and  q(x)  is  the  local  fixed 

B— 8 

charge  density  which  is  measured  by  a  CV  profiler.  To  measure 

the  trap  emission  rate  at  a  single  field,  two  trap  filling  pulses 
are  used.  These  are  shown  in  Figure  B-l.  The  first  pulse  is 
slightly  larger  than  the  second  and  it  brings  free  carriers  into 
traps  in  a  slightly  higher  field  region  than  the  second  pulse. 

Data  are  taken  in  such  a  way  as  to  subtract  the  emission  from  the 
traps  filled  by  the  two  pulses.  The  difference  is  due  to  the  traps 
filled  by  the  larger  pulse  and  not  filled  by  the  smaller  pulse. 
These  traps  are  in  a  region,  x,  which  can  be  determined  from  the 
profiler  data  for  x  versus  V  and  in  an  electric  field,  F(x),  which 
is  calculated  from  Equation  (B-l) .  These  effects  are  also  shown  in 
Figure  B-l.  By  using  various  pulse  heights,  the  position,  x,  can 


B  8G.L.  Miller,  IEEE  Trans.  Electron.  Devices,  ED-19,  1103  (1972). 
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be  moved  from  near  the  depletion  region  edge  to  near  the  junction 
and  fields  from  zero  to  ^  F  max  can  be  studied. 
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Figure  B-l.  Pulse  schematic  for  double  pulse  DLTS 
(a)  shows  the  applied  voltage  versus 
time,  (b)  shows  the  regions  filled  by 
the  two  pulses,  and  (c)  shows  the 
capacitance  signal.  The  DLTS  data  are 
the  difference  between  the  voltage 
taken  at  and  V2- 
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RESULTS 


A  DLTS  spectrum  for  Czochralski-grown  indium-doped  silicon  is 
shown  in  Figure  B-2 .  The  peak  directions  (all  upward)  indicate 
that  the  traps  are  hole  traps.  Four  traps  are  observed  and  they 
are  labeled  X,  In,  H^,  and  .  The  trap  concentrations  for  the 
sample  shown  are: 


In-X 
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Figure  B-2.  The  DLTS  spectrum  obtained  for 
indium-doped  silicon. 
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In  equilibrium  in  neutral  semiconductor  material  there  are  equal 
numbers  of  traps  thermally  emitting  trapped  carriers  as  there  are 
empty  traps  capturing  carriers.  This  equality  in  p-type  silicon 
can  be  used  to  relate  the  trap  emission  time  constant  to  its  cap¬ 
ture  coefficient,  giving 


(B-2) 


T 


e 


where 

Tg  =  Trap  emission  time  constant 

Erp  =  Trap  energy  relative  to  the  valance  band 

C  =  Capture  coefficient 

g  =  Ground  state  degeneracy 

m^  =  Hole  density  of  states  effective  mass 

h  =  Plank's  constant 

Using  the  relation 

a  =  C<v>  (B-3) 

the  average  thermal  velocity  <v>  =  (8KT/7rm*)^^  can  be  calculated 
and  Equation  (B-2)  can  be  put  in  the  form 
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(B-4) 


T2t 


e 


Et/KT 


2K“ 


2.3 

TT  h  g 


mh 


In  this  equation  the  density  of  states  and  velocity  effective 
masses  have  been  approximated  as  being  equal.  Both  the  In-X  and 
the  In  centers  show  a  strong  field  dependence.  The  double-pulse 
technique  was  used  to  measure  the  trap  emission  rate,  x,  as  a 
function  of  electric  field  and  temperature.  The  data  for  sub¬ 
stitutional  indium  are  shown  in  Figure  B-3,  where  it  is  plotted 
as  InT  x  versus  1/KT  to  fit  Equation  (B-4).  In  this  form  the 
data  should  be  in  a  straight  line,  with  the  slope  equaling  the 
effective  thermal  activation  energy. 


The  electric  fields  in  the  depletion  region  lower  the  trap  activa¬ 
tion  energies  and  increase  the  trap  emission  rate.  This  is  called 
the  Poole-Frenkel  effect  and  is  shown  in  Figure  B-4  for  a  coulomb- 
attractive  trap.  The  change  in  emission  rate  for  a  coulomb  trap 
in  an  electric  field  has  been  calculated  by  Hartke  as  being6-9 


Te(0) 

xe(F) 


1  + 


re/F 

[  KT 


1]  e  8/F/KT 


+  1/2 


(B-5) 


where 

x  (0)  =  Zero  field  emission  rate 
x  (F)  =  Emission  rate  in  a  field,  F 
6  =  (e3/Tre)1/2 

B-9 

J.L.  Hartke,  J.  Appl .  Phys .  39,  4871  (1968). 
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Figure 


V/m  to  2.5  x  106  V/m  are  shown  by  the  solid 
lines  in  the  lower  right.  The  dashed  lines 
are  corrections  of  these  data  to  zero  field. 
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NO  APPLIED  FIELD  F  *0  (a) 


Figure  B-4.  The  thermal  activation  energy  of  a  coulombic 
trap  in  zero  field  is  shown  in  (a).  The 
change  in  trap  energy  (AET),  the  lowered 
activation  energy  (Em),  and  the  thermally 
assisted  tunneling  eifect  are  shown  for  an 
applied  electric  field  in  (b). 

This  is  based  on  a  three-dimensional  model.  A  simpler  formula 

B—  in 

based  on  a  one-dimensional  model  is  often  used  instead  of  this. 

B/F 

t(0)  =  x(F)e  KT  (B-6) 


Another  effect  which  increases  the  emission  rate 


in  a  field  is 


thermally  assisted  tunneling.  This  is  also  shown 


in  Figure  B-4 . 


B  Frenkel,  Phy .  Rev .  54,  647  (1938). 
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A  one-dimensional  approximation  has  been  developed  to  include  this 

B— 11 

by  Vincent,  Chantre  and  Boise 


t(0)  =  t(F) 


AE  /KT  r E  /KT 
e  +  /  exp 

Jaet/kt 


Z-Z3/2 


I  i  h\3/2 

1 4  (2m*)2  ( KT ) 2  I 

| 3  qhF  /  x 


I  AEt|5/3 
I  ZKT  J 


where  E^,  is  the  trap  energy  and  AE^,  is  8/F. 


(B-7) 


We  have  modified  this  to  a  three-dimensional  form  by  defining  an 
effective  three-dimensional  barrier  change,  AE*,  by 


AE*/KT 


KT  | 

8/P 

2  | 

1  + 

KT 

-1 

e8/F/KT 

i 

+  1/2 


(  B-8  ) 


The  data  in  Figure  B-3  have  been  corrected  to  zero  field  values 
using  the  different  models:  the  one-dimensional  Poole-Frenkel 
model  [Equation  (B-6)],  the  three-dimensional  Poole-Frenkel  model 
[Equation  (B-5)],  the  one-dimensional  model,  including  tunneling 
[Equation  (B-7)],  and  the  three-dimensional  model,  including  tun¬ 
neling  [Equation  (B-8)].  Data  for  the  In-X-center  indium,  and  the 
two  deeper  traps  are  shown  in  Figure  B-5.  The  data  for  the  indium- 
X  center  and  indium  have  been  corrected  for  the  three-dimensional 
Poole-Frenkel  effect  and  thermally  assisted  tunneling  by  Equation 
(B-8).  The  activation  energies  measured  were  0.095  eV  for  In-X 
and  0.126  eV  for  indium.  The  correction  for  electric  field  effects 


B  11G.  Vincent, 
(1979) . 


A.  Chantre, 


and  D.  Bois,  J.  Appl .  Phys .  50,  5484 
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I/T  X  1000 


25  12.5  8.33  6.25  50 


o 

Figure  B-5.  The  T  x  data  are  shown  for  the  H^,  H2,  In,  and  In-X 
centersein  silicon  in  an  electric  field  of  -8.6  x  105 
V/m.  The  data  for  the  In  and  In-X  centers  have  been 
corrected  to  zero  field  by  the  three-dimensional 
coulomb  correction,  including  thermally  assisted 
tunneling.  These  are  the  dashed  lines. 


in  Equation  (B-8)  gave  corrected  energies  of  0.105  eV  for  In-X  and 
0.143  eV  for  indium.  The  experimental  data  for  the  deeper  levels 
gave  Ev  +  0.31  eV  and  Ev  +  0.45  eV  for  and  Hg,  respectively. 
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The  DLTS  techniques  can  be  used  to  determine  the  capture  coeffi¬ 
cient  as  a  function  of  temperature.  For  very  short  filling  pulses, 
the  free  carriers  are  moved  in  and  out  of  the  depletion  region  too 
quickly  to  be  trapped  and  no  DLTS  peak  is  observed.  For  very  long 
fill  pulses,  all  of  the  traps  saturated.  The  peak  amplitude  as  a 

function  of  filling  time  increases  exponentially  with  a  time  con- 
B-5 

stant  equal  to: 


T 


-1 

f 


Cp  =  a  <v>p 


( B-9 ) 


where  p  is  the  hole  free-carrier  concentration. 


Data  have  been  taken  to  determine  the  capture  coefficients  for 
the  traps  observed.  The  results  for  In  and  In-X  are  shown  in  Fig¬ 
ure  B-6 .  There  is  a  slight  temperature  dependence  for  indium  in 
this  temperature  region  and  very  little  temperature  dependence 
for  the  In-X  center.  The  data  can  be  fit  to  an  exponential  form 
with 


'In 


„  „  in-9  3.  +0.031  eV/KT 

=  7.6  x  10  cm  /sec  e 


and 


CT  =  7.7 
In-x 


ho-8  3.  +0.006  eV/KT 

x  10  cm  /sec  e 


(B-10 ) 


(B-ll) 


The  accuracies  of  the  activation  energies  are  'v  ±0.02  eV  for 
indium  and  +0.005  eV  for  indium  X.  The  free-carrier  thermal 
velocity  has  been  estimated  using  v  =  v/8KT/ Trm*  with  temperature- 
dependent  effective  mass  data  supplied  by  J .  Lang  of  the  Air 
Force  Materials  Laboratory.  This  allows  the  capture  cross  section 
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I/KT  (ev1) 


Figure  B-6.  Capture  coefficients  as  a  function  of 

temperature  for  indium  and  the  indium-X 
centers  in  silicon. 


to  be  calculated  as  a  =  C/<v>. 
coefficients,  and  capture  cross 
X  center  near  77  K  are  given  in 


Values  of  the  fill  times,  capture 
sections  for  indium  and  the  indium- 
Table  B-l . 
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TABLE  B-l .  TRAPPING  PARAMETERS  FOR 
In  AND  In-X 


Indium 

Temp . 

(K) 

Tf 

(ns) 

P 

,  -3. 

(cm  ) 

<v> 

(cm/sec) 

C 

,3.  . 

(cm  /sec) 

c 

,  2. 

(cm  ) 

96.7 

30 

1014 

1.08  x  107 

3.17  x  10“  7 

-14 

2.9  x  10 

77.4 

12.5 

1014 

9.84  x  106 

8  x  10-7 

-14 

8.13  x  10 

Indium- 

X 

77.4 

50 

1014 

9.84  x  106 

2  x  10-7 

2.03  x  10-14 

61.1 

38 

1014 

8.89  x  106 

2.6  x  10-7 

2.92  x  10-14 

The  capture  coefficients  for  the  two  deep  levels,  and  were 

estimated  by  fitting  Equation  (B-2).  This  assumes  that  there  is  no 
temperature  dependence  for  their  capture  coefficients  and  that 
field  effects  on  these  centers  are  small.  The  g  value  used  in 
Equation  (B-2)  was  taken  to  be  2,  which  is  what  is  expected  for 
a  nontetrahedral  defect  level.  The  values  obtained  are: 

CL.  =  2.0  x  10-^  cm^/sec  (B-12) 

H1 

and 

C„  =  1.2  x  10  **  cm^/sec  (B-13) 

n2 
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DISCUSSION 


The  DLTS  spectrum  shown  in  Figure  B-2  is  able  to  distinguish  four 
hole  traps  in  Czochralski-grown  indium-doped  silicon.  The  In-X 
center  is  very  strong  in  the  data  shown;  it  is  not  always  observed. 
Several  models  for  the  X-center  have  been  proposed.2-1,  B-3,  2-4 
Evidence  supporting  a  nearest-neighbor  carbon-indium  pair  model 
has  been  presented  by  this  group  and  Baron  et  al.2-3,  2-4 »  2-12 
The  identities  of  the  deeper  centers  and  Hg  are  not  known . 

Both  the  indium  and  the  indium-X  center  show  the  strong  electric 
field  effects  expected  for  coulomb  traps.  The  electric  field  in 
the  diode  depletion  region  lowers  the  effective  trap  activation 
energies  and  increases  their  emission  rates.  This  was  shown  in 
Figures  B-3  and  B-5,  where  a  2.5  x  10  V/m  electric  field  is  shown 
to  change  the  indium  emission  rate  by  one  to  two  orders  of  mag¬ 
nitude  and  lower  the  indium  activation  energy  by  0.035  eV.  The 
data  for  indium  shown  in  Figure  B-3  was  corrected  for  this  field 
effect  using  one-  and  three-dimensional  models.  There  is  a  signif¬ 
icant  difference  between  using  the  simple  one-dimensional  coulomb 
barrier  lowering  corrections  of  Equations  (B-6)  and  (B-7)  and  the 
more  accurate  three-dimensional  form  in  Equations  (B-5)  and  (B-8). 
The  emission  rates  for  the  two  corrections  differed  by  a  factor 
of  two  to  six  over  the  temperature  range  studied  and  the  activation 
energies  differed  by  ^  0.02  eV.  Because  of  this,  the  more  accurate 
three-dimensional  correction  is  used  for  the  In-X  and  In  data 
shown  in  Figure  B-5. 


B— 1 2 

C.E.  Jones,  D.  Schafer,  W.  Scott,  and  R.J.  Hager  in  the  report 
of  the  IRIS  Specialty  Group  on  Infrared  Detectors,  June  1979, 
Minneapolis,  MN. 
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The  data  reported  are  for  samples  with  electric  fields  in  the 
5  6 

10  V/m  to  10  V/m  range.  For  these  fields  and  for  the  trap  depths 
studied,  the  effect  of  thermally  assisted  tunneling  is  small.  The 

tunneling  correction  is  expected  to  be  more  important  at  higher 

,  B-12 

fields . 


The  capture  coefficients  for  indium  and  the  indium-X  center  were 
measured  by  the  trap  filling  method  of  Equation  (B-9).  This  tech¬ 
nique  involves  the  trap  filling  when  the  depletion  region  is  de¬ 
creased  and  it  occurs  in  field  free  conditions.  The  indium  trap 

filling  occurred  so  quickly  even  in  material  doped  to  as  low  as 
14  3 

10  B/cm  that  there  is  appreciable  scatter  to  the  data  shown 

in  Figure  B-6 .  The  indium  capture  coefficient  at  77K  is  given  as 

8  x  10  cm  /sec.  The  activation  energy  is  estimated  at  0.03  eV 

±0.02  eV  for  indium.  Blakemore  used  a  photoconductivity  technique 

-6  3 

to  determine  an  indium  capture  coefficient  of  ^  1  x  10  cm  /sec 

B—  1 3 

at  77K  with  an  activation  energy  of  ^  0.05  eV.  Forbes  used  a 

MOSFET  technique  similar  to  DLTS  and  measured  a  capture  coefficient 

r*  q 

of  0.4  to  1.7  x  10-  cm  /sec  at  77K  with  activation  energy  of  0.05 
R— 14 

eV.  Forbes  used  bulk  values  of  p  for  his  calculations.  Fields 

in  the  depletion  region  empty  traps  in  this  region,  making  the 
fixed  charge  here  different  from  that  in  the  bulk  material .  When 
the  filling  pulse  is  applied,  carriers  move  in  to  neutralize  the 
local  fixed  charge,  making  p  in  Equation  (B-9)  equal  to  the  de¬ 
pletion  region  fixed  charge  not,  as  is  often  done,  equal  to  the 


B  13J.S.  Blakemore,  Can.  J.  Phys .  34,  (1956). 

B-14L.  Forbes  in  "Characterization  and  Analysis  of  Indium-Doped 
Silicon  Extrinsic  Detector  Material,"  (J.  Backus,  Hughes  Res. 
Lab.,  T.  McGill,  C.I.T.;  and  L.  Forbes,  U.  of  C.  Davis) 
Interim  Technical  Report  1,  Dec.  1977,  DARPA  Contract  DAAK70- 
77-C-0082 ,  pg.  48-51. 


65 


free-carrier  concentration  in  the  neutral  bulk  semiconductor.  The 
value  of  p  needed  for  Equation  (B-4)  can  be  measured  on  a  CV 
profiler  which  measures  the  fixed  charge  in  the  depletion  region. 
Using  a  depletion  region  value  of  p  lowers  Forbes  C  value  and  also 
lowers  his  activation  energy  which  would  bring  them  closer  to  the 
values  determined  in  this  study. 

The  average  capture  coefficient  value  measured  by  DLTS  and  reported 
herein  is  very  close  to  the  values  reported  by  Blakemore  and  Forbes. 
The  indium  capture  coefficient  activation  energy  of  0.03  ±0.02  eV 
has  a  large  error  range,  but  within  the  error  range  it  is  the  same 
as  that  reported  by  the  other  two  studies.  This  gives  some  con¬ 
fidence  in  the  DLTS  data  for  the  indium-X  center  capture  coeffici¬ 
ent  which  has  much  less  experimental  error. 

The  magnitude  of  the  indium  and  indium-X  capture  coefficients  of 
-6  -7  3 

10  to  10  cm  /sec  are  both  large  and  in  the  range  expected  for 

coulomb  traps.  The  capture  coefficients  estimated  for  H1  and  Hn 

-8  -9  12 
were  smaller  in  the  10  to  10  range,  which  suggests  that  they 

may  be  neutral  traps.  The  effect  electric  fields  have  on  neutral 

traps  is  much  less  than  that  on  coulombic  centers.  The  capture 

coefficients  calculated  from  Equation  (B-2)  for  the  H1  and  H0 

traps  neglected  any  field  dependence. 

Trap  energies  measured  by  DLTS  can  be  corrected  for  electric  field 
effects  and  the  temperature  dependence  of  the  capture  cross  sec¬ 
tion.  Data  to  make  both  of  these  corrections  have  been  measured 
for  the  In  and  In-X  centers,  as  follows: 
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Center 

edlts 

+EField 

-E 

a 

ETrap 

Optical 

Trap 

In 

0.126  eV 

+  0.017  - 

0.03 

=  0.113  eV 

0.156  eV 

In-X 

0.095  eV 

+  0.010  - 

0.006 

=  0.099 

0.1128  eV 

In  these  cases  we  are  adding  three  terms  each  with  an  error  of 
^  0.01  eV  and  coming  out  low  by  0.04  eV  to  0.014  eV  for  In  and 
In-X.  Again,  the  In-X  data  had  less  experimental  error  for  E 
and  the  corrected  DLTS  trap  energy  is  closer  to  the  optical  value. 
The  field  correction  used  also  may  be  low.  Using  the  one-dimensional 
formulas  Equations  (B-6)  and  (B-7)  would  give  higher  field  correc¬ 
tions  than  the  three-dimensional  formulas.  Better  correlation 
between  theory  and  experiment  does  not  change  the  fact  that  the 
one-dimensional  formulas  were  based  on  a  less-accurate  physical 
model,  however.  There  is  a  large  experiment  error  in  the  expo¬ 
nential  dependence  of  the  capture  coefficients  and  these  terms 
could  be  smaller,  bringing  the  optical  and  DLTS  energies  closer. 

CONCLUSION 

The  DLTS  technique  has  identified  four  traps  in  indium-doped 

silicon.  Two  of  these  are  the  indium-X  center  and  substitutional 

indium  and  two  are  unidentified  deeper  centers  at  Ev  +  0.31  eV  and 

E  +  0.45  eV.  Both  indium  and  the  indium-X  center  show  the  large 
v 

electric  field  effects  expected  for  coulomb  traps.  The  electric 
field  increases  trap  emission  rates  and  decreases  their  activation 
energies.  The  emission  rates  and  activation  energies  have  been 
measured  for  a  range  of  electric  fields.  Theoretical  formulas  have 
been  able  to  correct  these  data  to  single-zero-field  curves  for 
indium  and  the  indium-X  center.  The  resulting  zero-field  trap 
activation  energies  are  lower  than  the  optically  measured  values, 
however . 
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The  low-temperature  capture  coefficients  for  the  four  hole  traps 
have  been  measured . Indium  is  the  only  one  of  these  which  has  been 
previously  reported  in  the  literature  and  the  DLTS  values  measured 
agree  with  the  reported  value  for  this  trap.  The  indium  hole  trap¬ 
ping  occurs  so  quickly  that  there  is  a  large  error  in  the  DLTS 
determination  of  the  temperature  dependence  of  the  capture  co¬ 
efficient  and  the  DLTS  corrected  trap  energy  which  uses  these  data. 
The  data  for  the  indium-X  center  have  less  experimental  error  but 
still  yield  a  trap  energy  which  is  slightly  lower  than  optical  or 
Hall  values. 

The  DLTS  data  provide  the  trapping  parameters  which  can  be  used 
with  optically  measured  energy  levels  to  mathematically  model  trap¬ 
ping  effects  in  indium-doped  silicon. 
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APPENDIX  C 


TRIGONAL  STRESS  SYMMETRY 
OF  THE  Al-X  CENTER  IN  SILICON* 

ABSTRACT 

Electrically  active  complexes  called  X-centers  have  been  identified 
associated  with  the  group  IIIA  acceptors  in  silicon.  The  complexes 
are  acceptors  with  an  energy  level  ^  80%  that  of  the  correspond¬ 
ing  group  IIIA  substitutional  element.  To  provide  the  information 
needed  to  help  identify  these  defects,  the  optical  absorption 
spectra  for  the  aluminum-X  center  have  been  measured  under  condi¬ 
tions  of  uniaxial  stress.  The  data  are  interpreted  in  terms  of  a 
defect  having  trigonal  [ill]  symmetry.  Defect  models  which  have 
been  suggested  involving  single  substitutional  atoms  or  distance 
pairs  do  not  fit  the  stress  data.  On  the  other  hand,  nearest- 
neighbor  substitutional  pairs  or  a  substitutional  atom  next  to  a 
bond-centered  or  tetrahedral  interstitial  fit  the  data. 

INTRODUCTION 

Recent  investigations  of  silicon  doped  with  group  IIIA  acceptors 
(B,  A1 ,  Ga,  and  In)  have  revealed  several  shallow  acceptor  levels 


*The  authors  of  this  paper  are  David  E.  Schafer,  Colin  E.  Jones, 
and  J.E.  Sjerven,  Honeywell  Corporate  Material  Sciences  Center, 
Bloomington,  MN.  The  paper  was  submitted  to  J.  Appl .  Phys .  for 
publication . 
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referred  to  as  "X"-levels.C  1  C  4  Each  X-level  appears  to  be  a 
complex  involving  a  particular  group  IIIA  impurity  but  with  an 
ionization  energy  20  to  25%  lower  than  the  corresponding  substitu¬ 
tional  group  IIIA  element  in  silicon.  The  structure  of  the  X- 
center  has  not  been  firmly  established.  However,  several  models 
for  X-centers  have  been  suggested,  including  centers  involving 
carbon-indium  pairs,  ’  indium  interstitial  or  vacancy  com- 

plexes,  and  indium-aluminum  pairs.  In  the  case  of  Si(In), 
the  X-defect  is  of  current  technological  interest  in  connection 
with  the  use  of  this  material  for  a  3-5  micron  photoconductor, 
since  the  shallower  In-X  level  contributes  excess  thermal  gen- - 
erat ion-recombination  noise  and  requires  lower  temperatures  to 
achieve  photon  noise-limited  operation. 

Toward  an  eventual  identification  of  the  X-defect(s),  we  have 
investigated  the  behavior  of  the  Al-X  optical  absorption  spectrum 
under  uniaxial  stress.  This  has  provided  information  about  the 
symmetry  and  orientation  of  the  Al-X  defect.  Although  the  X-defect 
in  Si(Al)  is  not  of  as  great  current  technological  importance  as 
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that  in  Si (In),  it  is  the  most  amenable  of  the  group  IIIA  X- 

centers  to  an  optical  study  of  this  kind,  having  a  relatively 

large  optical  cross  section  and  having  optical  transitions  easily 

separable  from  other  nearby  spectral  features.  Based  on  consistent 

C_  2  c  _  3  C— 4 

trends  in  the  ionization  energies  *  and  annealing  behavior 

of  the  X  defects  in  B-,  A1-,  Ga- ,  and  In-doped  silicon,  it  seems 

likely  at  present  that  the  various  X-defects  in  these  cases  have 

related  origins  and  that  the  structural  information  given  here  for 

the  Al-X  defect  will  be  applicable  to  the  other  X-defects  as  well. 

Experimental  Procedure 

1 6 

Oriented  samples  of  Czochralski-grown  silicon  with  3  to  5  x  10 
Al/cm“3  were  cut  with  dimensions  0.6  x  0.45  x  0.13  inch,  the  trans¬ 
mission  faces  angle-lapped  at  about  0.5°  to  prevent  interference 
effects,  and  polished.  They  were  mounted  in  a  copper  sample  holder 
and  compressed  along  the  [ill]  ,  [100]  or  [110]  axis  by  differen¬ 

tial  thermal  contraction  relative  to  copper  upon  cooling  to  10K, 
the  temperature  of  the  measurements .  Some  samples  were  mounted 
conventionally  for  stress-free  data.  Infrared  absorption  measure¬ 
ments  were  made  using  a  Digilab  FTS-14  Fourier  transform  spectrom¬ 
eter  at  a  resolution  of  2  cm"^.  In  some  of  the  measurements,  light 
polarized  either  parallel  or  perpendicular  to  the  stress  direction 
with  a  wire  grid  polarizer  was  used. 

Experimental  Results 

The  optical  absorption  spectrum  for  unstressed  aluminum-doped 
silicon  in  the  200-cm-1  to  450  cm-1  region  is  shown  in  Figure 
C-l.  Data  are  shown  for  sample  temperatures  of  9K  and  50K.  The 
peaks  at  245,  278,  and  320  cm-1  are  due  to  substitutional  boron 
which  is  present  in  this  sample  as  an  impurity.  The  peaks  at 
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ABSORBANCE 


339,  372,  404,  and  414  cm  ^  are  due  to  the  aluminum-X  center. 

The  general  rise  in  the  spectra  toward  larger  wave  numbers  is  due 
to  the  strong  substitutional  aluminum  absorption  which  has  its 
first  peak  at  440  cm-'*'. 


300 


400 


FREQUENCY  (cm1) 


Figure  C-l.  The  absorption  spectrum  of  Czochralski- 
grown  aluminum-doped  silicon,  showing 
lines  due  to  the  Al-X  center  (XI,  X2, 
X3,  X4)  and  due  to  residual  boron  (Bl, 
B2,  B4).  Spectral  resolution  is  2  cm-1. 
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The  absorption  spectra  with  [100]  and  [ill]  uniaxial  stress  applied 
are  shown  in  Figures  C-2  and  C-3,  respectively.  These  data  are 
taken  at  10K  with  light  polarized  parallel  and  perpendicular  to 
the  stress  directions. 

For  stress  applied  in  the  [110]  direction  there  are  two  inequiva¬ 
lent  directions  for  the  optical  beam.  Data  for  [HO]  uniaxial 
stress  with  the  optical  beam  in  the  [110]  and  [00l]  directions 
are  shown  in  Figure  C-4  and  C-5,  respectively. 

Theory 

The  stress  splittings  and  degeneracies  of  hydrogenic  acceptor 
levels  in  silicon  as  predicted  on  the  basis  of  symmetry  properties 
are  summarized  in  Figure  C-6 .  These  predictions  are  the  basis  for 
the  analysis  of  the  stress  spectra  presented  below.  The  simplest 
picture  of  an  acceptor  in  silicon,  that  of  a  positive  hole  at¬ 
tracted  to  a  fixed  negative  ion  via  a  spherically  symmetric 
potential,  is  represented  at  the  left  of  the  figure.  It  consists 

of  a  ls-like  ground  state  with  allowed  optical  transitions  to 

C—8 

three  2p-like  excited  states.  Since  shallow  acceptor  levels 
generally  arise  from  the  four-fold  degenerate  J3/2  valence 
band  maximum,  the  Is  and  2p  hydrogenic  states  at  the  left  of 
Figure  C-6  are  assigned  net  degeneracies  of  4  and  12,  respectively. 
The  level  diagrams  in  the  remainder  of  the  figure  represent  the 
effects  of  various  deviations  from  spherical  symmetry,  due  to 
local  structure  or  applied  stress.  In  this  figure,  these  effects 
are  represented  as  small  perturbations  on  the  ls-2p  hydrogenic 
level  structure. 


C—8 

W.  Kohn,  in  Solid  State  Physics,  Vol.  5,  p.  257,  F.  Setz  and 
D.  Turnbull,  Eds.  (Academic  Press,  New  York,  1957). 
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Figure  C-2.  The  absorption  spectrum  of  the  material 
of  Figure  C-l  under  [100]  uniaxial  com¬ 
pression  near  liquid-helium  temperature 
The  inferred  Al-X  level  splittings  are 
also  shown.  The  spectral  resolution  is 
2  cm-l. 
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Figure  C-3.  The  absorption  spectrum  of  the  material  of 

Figure  C-l  under  [ill]  uniaxial  compression. 
The  inferred  Al-X-level  splittings  are  also 
shown,  with  two  configurationally  inequiva¬ 
lent  subsets  of  the  centers  represented  in 
two  separate-level  diagrams . 
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Figure  C-4.  The  absorption  spectrum  of  Al-doped  silicon 
under  [HOl  uniaxial  compression  with 
incident  light  along  [110]  .  The  inferred 
level  splittings  are  also  shown,  with  two 
inequivalent  subsets  of  the  centers  repre¬ 
sented  in  two  separate-level  diagrams. 
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Figure  C-5.  The  absorption  spectrum  of  Al-doped  silicon 
under  [lio]  uniaxial  compression  with  the 
incident  light  along  [001] .  The  inferred 
level  splittings  are  also  shown,  with  two 
inequivalent  subsets  of  the  centers  repre¬ 
sented  in  two  separate-level  diagrams . 
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Figure  C  6.  Splitting  patterns  of  Is  and  2p  hydrogenic  acceptor 
levels  under  the  influence  of  local  defect  struc¬ 
ture  and  applied  stress,  as  predicted  on  the  basis 
of  symmetry  changes.  The  upper  half  of  the  figure 
treats  the  case  of  a  substitutional  acceptor 
(tetrahedral  symmetry),  while  the  lower  half  depicts 
the  case  of  a  trigonal  structure  point  group  C3V, 
with  the  threeef old  axis  oriented  along  any  of  the 
body  diagonals  of  the  cubic  host  structure. 


The  top  half  of  Figure  C-6  treats  the  usual  case  for  group  IIIA 
acceptors,  that  of  a  substitutional  impurity  in  a  tetrahedral 
environment.  As  indicated  in  the  top  center  of  Figure  C-6,  the 
local  tetrahedral  symmetry  is  expected  to  remove  the  12-fold 
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degeneracy  of  the  2p  excited  state,  splitting  it  into  four  compo¬ 
nents  of  degeneracy:  4,  4,  2,  and  2.  Such  splitting  and  degen¬ 
eracies  have  been  observed  for  the  most  part  in  investigations 

C-9 

of  the  absorption  spectra  of  the  group  IIIA  impurities,  with 

the  possible  exception  of  excited  state  4,  for  which  fourfold 

C-9 

degeneracy  has  been  suggested  rather  than  the  value  of  two 
predicted  by  the  effective-mass  theory.  The  energy  ordering  of 
the  fourfold  and  twofold  degenerate  states  in  Figure  C-6  has  been 
taken  from  the  results  for  the  substitutional  acceptors.  The 
additional  symmetry-related  effects  of  [100]  and  [ill]  applied 
stress  are  illustrated  in  the  upper  right  portion  of  Figure  C-6 
and  are  qualitatively  the  same  in  each  case,  resulting  in  six  2p 
levels  and  two  Is  levels  all  of  degeneracy  two.  This  is  the  maxi¬ 
mum  degree  of  splitting  which  these  levels  may  undergo,  since 
these  levels  have  only  the  spin  degeneracy  of  two  left.  Not  all 
of  the  possible  transitions  are  observed  with  the  same  intensity 
due  to  differing  selection  rules  for  the  various  transitions 
and  polarizations. 

Complex  defects,  as  the  X-centers  may  be,  have  a  symmetry  lower 
than  tetrahedral.  We  describe  the  case  of  a  trigonal  [ill]  de¬ 
fect  in  some  detail  below  and  then  sketch  the  effects  expected 
for  other  defect  symmetries.  A  trigonal  defect  has  a  point  group 
(Cgv  or  Dgd)  which  is  a  subgroup  of  the  tetrahedral  group.  The 
ground  state  is  split  by  the  nontetrahedral  crystal  field  into 
two  ground-state  spin— doublets .  The  excited  states  are  split  into 
six  components  as  shown  in  the  lower  left  of  Figure  C-6.  The  ex¬ 
cited  states  are  spatially  diffuse  and  have  a  node  at  the  defect, 


C-9A.  Onton,  P.  Fisher,  and  A.K.  Ramdas ,  Phys .  Rev.  163,  1967. 
(See  Figs.  7  and  25  Al-X  absorption  spectra. ) 
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so  the  excited  state  splittings  are  not  expected  to  be  very 
large.  The  ground-state  components  are  more  strongly  localized 
at  the  defect  and  they  could  be  appreciably  split . 

The  lower  right  hand  side  of  Figure  C-6  shows  the  splitting  ex¬ 
pected  under  applied  uniaxial  stress.  The  trigonal  [ill]  complex 
has  an  added  stress  effect  over  that  of  a  substitutional  tetrahedral 
defect.  A  trigonal  defect  can  be  oriented  along  any  of  the  four 
cube  diagonals.  The  crystal  will  contain  equal  numbers  of 
defects  in  each  of  these  orientations  if  the  stress  is  applied 
at  a  low  enough  temperature  so  that  the  defect  cannot  reorient. 

If  there  are  different  angles  between  the  four  orientations  and 
the  applied  stress,  then  the  defects  in  the  different  orientations 
will  be  stressed  different  amounts.  Optical  data  will  show  absorp¬ 
tion  from  defects  in  each  of  the  different  orientations,  giving 
rise  to  a  superposition  of  several  sets  of  stress  spectra.  This 
is  called  a  configurational  effect.  The  intensities  of  the  lines 
from  a  particular  orientation  will  not  necessarily  be  in  the  same 
proportion  to  those  of  the  other  orientations,  however,  since  the 
intensity  of  each  transition  changes  with  viewing  angle. 

The  configurational  effect  is  expected  to  occur  predominantly  in 
the  ground  state.  Again,  because  the  excited  states  are  spread 
out  and  have  a  node  at  the  defect,  the  stress  splitting  of  the 
excited  states  is  expected  to  arise  mainly  from  the  distortion  of 
the  bulk  silicon  lattice  and  is  not  expected  to  be  strongly  af¬ 
fected  by  the  local  atomic  symmetry  of  the  defect.  Consequently, 
the  stress-induced  changes  in  these  excited  states  (bottom  half 
of  Figure  C-6)  are  represented  in  the  figure  as  being  the  same  as 
those  of  the  tetrahedral  defect  (upper  half  of  Figure  C-6). 
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Under  [lOO]  stress,  the  different  orientations  for  the  trigonal 
[ill]  defect  all  make  equal  angles  with  the  stress  and  will  be 
affected  in  the  same  way,  producing  a  single  stress  pattern  in  the 
optical  data.  For  a  [ill]  stress,  one  of  the  defect  orientations 
is  along  the  stress  direction  and  three  are  at  an  angle  to  it, 
which  would  be  expected  to  give  a  superposition  of  two  stress  pat¬ 
terns  . 

Application  of  [lio]  uniaxial  stress  also  creates  two  distinct 
subsets  of  [ill] -oriented  centers,  one  with  the  trigonal  axis 
perpendicular  to  the  stress  (two  possible  orientations)  and  the 
other  with  the  trigonal  axis  not  perpendicular  to  the  stress  (two 
more  orientations).  The  resulting  two  level  diagrams  will  be 
similar  to  those  shown  for  the  [ill]  stress  case. 

Similar  sets  of  level  diagrams  can  be  drawn  illustrating  the  stress 
effects  expected  for  rhombic  defects,  which  includes  the  [110] - 
and  [lOO]  -directed  centers,  and  monoclinic  defects.  The  patterns 
of  splittings  and  degeneracies  in  these  diagrams  are  similar  to 
those  for  the  trigonal  defect  illustrated  here.  The  unique  feature 
of  the  trigonal  [111]  defects  is  the  lack  of  configurational  ef¬ 
fects  for  a  [lOO]  stress. 

DISCUSSION 

Figure  C— 1  shows  the  Al-X  absorption  spectrum  in  unstressed 
silicon  at  two  temperatures.  Three  strong  lines  and  one  weak  line 
(labeled  XI  through  X4)  are  observed,  as  well  as  three  weak  lines 
(labeled  Bl,  B2,  and  B4)  due  to  residual  boron.  As  noted  else¬ 
where,*"-2  the  four  Al-X  lines  are  quite  similar  in  their  spacing 
and  relative  intensity  to  corresponding  lines  in  the  optical 
spectra  of  other  shallow  acceptors,  suggesting  that  an  effective- 
mass  picture  is  appropriate  for  the  Al-X-excited  states .  As 
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mentioned  above,  this  is  not  surprising,  since  these  states  are 
not  expected  to  interact  strongly  with  the  structural  details  at 
the  defect. 

The  local  crystal  field  of  any  nontetrahedral  defect  would  split 
the  ground  state  into  two  spin-doublet  levels.  The  possible 
presence  of  ground-state  splitting  due  to  nontetrahedral  defect 
symmetry  was  investigated  through  the  second  spectrum  of  Figure 
C— 1 ,  taken  at  50K.  At  this  temperature,  any  additional  ground- 
state  components  less  than  about  10  meV  above  the  observed  ground 
state  would  have  been  thermally  populated  and  revealed  through 
additional  optical  transitions.  The  absence  of  any  additional  lines 
implies  that  any  split-off  ground-state  component,  if  it  exists, 
is  either  less  than  2  meV  away  from  the  lower  ground  state  (i.e., 
unresolvable )  or  more  than  10  meV  above  it. 

The  apparent  effective-mass  nature  of  the  Al-X-excited  states  was 
used  in  the  analysis  of  the  Al-X  spectrum  under  applied  stress. 

It  was  assumed  that  there  would  be  a  close  correspondence  between 
the  stress  behavior  of  these  excited  states  and  that  of  the  (ef¬ 
fective-mass-like)  excited  states  of  other  shallow  acceptors.  Thus, 
the  stress-induced  changes  in  the  nearby  boron  lines  in  Figure  C-l 
were  used  as  a  guide  in  attributing  some  of  the  observed  Al-X  line 
splittings  to  the  Al-X  excited  states.  Also,  configurational  ef¬ 
fects  (i.e.,  multiple  orientations  of  A  nontetrahedral  defect 
relative  to  the  stress  direction)  were  neglected  for  the  excited 
states,  again  because  of  their  apparent  weak  interaction  with  struc¬ 
tural  details  near  the  defect. 

The  Al-X  spectrum  under  [100]  stress,  shown  in  Figure  C-2,  is  the 
simplest  of  the  three  stress  cases  to  discuss,  since  ground-state 
splitting  does  not  appear  to  play  a  role;  the  observed  spectral 
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changes  can  be  accounted  for  by  excited-state  splitting  only .  In 
this  case  we  estimate  the  expected  Al-X  excited-state  splitting 
from  the  behavior  of  nearby  boron  lines  Bl  and  B4 .  Previous  inves¬ 
tigations  have  shown  that  the  upward  shifts  in  these  lines  are  due 

C— 9 

to  splitting  in  the  boron  ground  state.  Split-off  downshifted 
partners  to  lines  Bl  and  B4  observed  in  the  earlier  work  are 
not  present  in  the  spectrum  of  Figure  C-2  because  of  depopulation 
of  the  upshifted  split-off  ground-state  component  associated  with 
those  transitions.  Published  values  of  excited-state  to  ground- 
state  splitting  ratios  for  boron  acceptors  give  splittings  of 
0,  14,  and  0  cm-'1'  for  boron-excited  states  1,  2,  and  4,  respec¬ 
tively  . 

Returning  to  the  Al— X  absorption  lines  we  find  line  splittings  of 
0,  13,  and  0  cm-1  for  lines  1,  2,  and  4,  respectively  which  is  in 
good  correspondence  with  the  estimated  boron-excited-state  split¬ 
tings.  It  is  thus  natural  to  attribute  all  of  the  observed  Al-X 
line  splittings  to  the  excited  states.  No  lines  are  observed 
which  can  be  attributed  to  ground-state  splitting.  The  inferred 
level  structure  in  the  stressed  condition  is  shown  in  the  upper 
part  of  Figure  C-2. 

Al-X  centers  which  have  different  orientations  relative  to  the 
applied  stress  would  be  expected  to  have  slightly  different 
ground-state  energies.  Any  ground-state  energy  differences  due 
to  these  configurational  effects  should  produce  extra  lines  in 
the  spectrum.  Depopulation  of  the  higher— energy— defect  orienta¬ 
tions  should  not  be  possible  in  these  samples.  They  are  not  highly 
compensated  and  the  applied  light  does  not  empty  any  appreciable 
fraction  of  the  centers.  There  are  thus  not  enough  ionized  cen¬ 
ters  available  to  depopulate  the  higher-energy  configurations. 
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Since  no  extra  Al-X  lines  are  observed  for  the  [100]  stress,  there 
appears  to  be  no  configurational  effect  for  this  stress  direction. 

The  [111]  stress  behavior  (Figure  C-3)  shows  a  higher  degree  of 
line  splitting  than  the  [100]  stress  case.  The  same  procedure 
described  for  [100]  stress  was  followed  in  order  to  identify  the 
Al-X  excited-state  splittings.  However,  in  order  to  account  for  all 
of  the  observed  transitions,  the  assumption  of  two  ground-state 
components  separated  by  about  6.5  cm-1  was  required. 

The  polarization  dependence  of  some  of  the  absorption  lines  in 
Figure  C-3  reveal  differences  from  the  stress  behavior  of  tetra¬ 
hedral  substitutional  acceptors.  Specifically,  the  appearance  of 
three  distinct  components  of  line  2  in  the  E | | F  polarization  is 
inconsistent  with  the  polarization  selection  rules  for  a  tetra¬ 
hedral,  defect;  the  two  stress-induced  ground-state  components  in 
the  tetrahedral  case  could  undergo  at  most  only  one  E | | F  transition 
each  to  any  component  of  excited  state  2,  permitting  at  most  two 
components  of  line  2.  Defect  structures  of  symmetry  lower  than 
tetrahedral,  on  the  other  hand,  involve  fewer  selection  rules  for 
optical  transitions  and  would  not  be  inconsistent  with  the  spectrum 
of  Figure  C-3.  Again,  the  inferred  arrangement  of  levels  in  the 
stressed  condition  is  shown  in  the  figure. 

Two  aluminum-X-center  spectra  under  [HO]  uniaxial  stress  are 
shown  in  Figures  C-4  and  C-5.  Analysis  of  observable  changes  in 
lines  2  and  4  as  before  (changes  in  line  1  are  partly  obscured 
by  boron  line  4)  indicate  two  ground-state  components. 

Considering  now  the  results  of  the  three  stress  spectra  collec¬ 
tively,  we  conclude  that  the  Al-X  defect  is  most  likely  to  be  of 
trigonal  symmetry,  with  four  possible  orientations  of  the  tri¬ 
gonal  axis  along  the  body  diagonals  of  the  cubic  host  structure. 
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The  number  of  ground-state  levels  occurring  for  each  direction 
of  applied  stress  [one  for  [100]  stress  and  two  for  [ill]  and 
[110]  stress]  corresponds  to  the  number  of  inequivalent  orienta¬ 
tions  expected  for  a  defect  of  such  symmetry. 

Of  the  structures  which  have  been  proposed  to  date  for  X-centers, 
those  having  trigonal  symmetry  and  which  are  therefore  consistent 
with  our  conclusions  here  are  illustrated  in  Figure  C-7 .  These 
include  nearest-neighbor  pairs  of  substitutional  impurities,  a  sub 
stitutional-bond-centered  interstitial  pair  or  a  substitutional- 
tetrahedral  interstitial  pair.  Distant  pairs  of  impurity  atoms  or 
single  substitutional  impurities,  on  the  other  hand,  are  not  con- 


Figure  C-7.  Three  defect  structures  having  the  trigonal, 
[ill] -oriented  symmetry  discussed  in  the 
text:  (a)  nearest-neighbor  impurity  pair; 

(b)  bond-centered  interstitial  impurity; 
and  (c)  a  substitutional-tetrahedral  inter¬ 
stitial  pair. 
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sistent  with  the  trigonal  symmetry  identified  here.  The  association 
of  the  (presumably  analogous)  In-X  defect  with  the  presence  of 
carbon  impurities  tends  to  rule  out  interstitial  aluminum  as  a 
possible  Al-X  structure. 
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APPENDIX  D 

CARBON- ACCEPTOR  PAIR  CENTERS 
(X-CENTERS)  IN  SILICON* 


An  electrically  active  complex  called  the  X-center  has  long  been 
known  in  indium-doped  silicon.  Similar  defects  have  now  been  asso¬ 
ciated  with  the  other  group  IIIA  elements  in  silicon.  These  X- 
centers  are  acceptors  with  energies  approximately  80%  of  the 

corresponding  group  IIIA  substitutional  atom's  energy  and  with 

-2  -5 

concentrations  10  to  10  of  the  corresponding  group  IIIA  con¬ 
centration.  Doping  studies  show  that  the  X-center  concentrations 

depend  on  carbon  and  fit  the  reaction  A  +  C  -*•  X,  where  A  is  the 

s  s  s 

substitutional  group  IIIA  acceptor,  C  is  the  substitutional 

s 

carbon,  and  K  is  the  equilibrium  constant.  The  equilibrium  con¬ 
stants  as  a  function  of  temperature  and  the  complex  binding  energ¬ 
ies  have  been  determined  for  all  of  the  group  IIIA  X-centers 
except  thallium.  The  binding  energy  for  the  complex  increases 
from  0.03  eV  for  the  boron-X  center  to  0.70  eV  for  the  indium-X 
center.  The  complexes  do  not  anneal  out  up  to  the  melting  point 
of  silicon  and  they  can  be  cyclically  annealed.  Stress  studies  on 
the  aluminum-X  center  show  that  the  complex  has  [ill]  trigonal 
orientational  symmetry.  The  dependence  on  carbon,  the  stress  sym¬ 
metry,  and  the  high-temperature  stability  of  these  centers  all 
suggest  a  model  for  these  complexes  of  a  nearest-neighbor  sub¬ 
stitutional  carbon-acceptor  pair.  Mathematical  modeling  using 
extended  Huckel  calculations  confirms  that  this  defect  should 


*Work  done  under  DARPA  contract  DAAK70-77-C-0194 .  The  authors  of 
this  paper  are  Colin  E.  Jones,  David  E.  Schafer,  Walter  M.  Scott, 
and  Robert  J.  Hager,  Honeywell  Corporate  Material  Sciences  Cen¬ 
ter,  Bloomington,  MN.  The  paper  was  submitted  to  J.  Appl .  Phys . 
for  publication. 
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have  a  single  energy  level  slightly  shallower  than  the  substitu¬ 
tional  acceptor. 

INTRODUCTION 


Carbon  and  oxygen  are  the  major  residual  impurities  in  silicon. 

Carbon  is  often  detected  at  concentrations  of  1017  c/cm2  in 

16  3 

Czocharlski  and  10  C/cm  in  float-zone  silicon.  Most  of  this 
carbon  is  in  isolated  substitutional  sites  and  is  electically  in¬ 
active.  This  paper  discusses  the  properites  of  a  class  of  elec¬ 
trically  active  defects  which  we  believe  to  be  carbon-acceptor 
complexes . 

Historically ,  a  shallow  acceptor  called  the  X— center  has  been  seen 

in  Hall  and  optical  absorption  dataD-2  in  indium-doped  silicon. 

This  center  is  known  to  involve  indium.  It  is  electrically  active, 

having  an  acceptor  level  at  +  0.11  eV,  and  it  occurs  in  pulled 

silicon  at  concentrations  0.1  to  0.5%  of  the  indium  concentrations. 

Silicon  crystals  doped  with  indium  or  other  group  IIIA  elements 

are  being  studied  as  infrared  photodetector  materials.  In  these 

devices,  even  at  the  low  concentration  levels  seen,  1012  x/cm2 
15  3 

to  10  X/ cm  ,  the  X-center  is  a  source  of  thermal  noise  in  the 

photodetectors.  To  eliminate  this  noise,  the  detectors  have  to  be 
cooled  to  freeze  out  the  X— centers.  This  involves  operating  the 
detectors  at  a  lower  temperature  than  would  be  needed  if  there 
were  no  X-levels  present.  In  some  applications,  this  extra  cooling 
is  a  severe  handicap. 


D  1R. Baron,  M.H.  Young  J.K.  Neeland,  and  O.J.  Marsh,  Appl.  Phvs 
Lett.  30,  594  (1977).  — — 

D-2 

M.W.  Scott,  Appl.  Phys.  Lett.  32,  540  (1978). 
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Several  models  have  been  suggested  for  the  indium-X  center.  The 

center  has  been  found  to  depend  on  carbon,  and  substitutional 

D-3  D-4 

indium-carbon  pair  model  has  been  suggested.  ’  The  energy 

level  of  the  defect,  between  indium  and  aluminum,  has  also  sug- 

D-5 

gested  In-Al  distant  pairs.  The  fact  that  the  X-center  concen¬ 
tration  increases  with  irradiation  has  suggested  that  the  center 
may  be  an  indium-interstitial  or  vacancy-related  defect.  Finally, 

EPR  data  have  shown  transition  metal  ions  and  acceptors  form 

D— 7  D— 8 

complexes.  This  has  suggested  Fe^-Ing  pairs  for  the  X-Center. 

This  paper  has  two  major  aims.  The  first  is  to  show  that  there  is 
an  X-center  for  each  of  the  group  IIIA  elements  in  silicon.  The 
second  is  to  present  the  optical,  thermal,  and  chemical  evidence 
we  have  obtained  which  supports  the  identification  of  these  de¬ 
fects  as  being  nearest-neighbor  substitutaional  carbon-acceptor 
pairs . 


D  3R.  Baron,  J.P.  Bankus ,  S.D.  Allen,  T.C.  McCill ,  M.H.  Young, 

H.  Kumura,  H.V.  Winston,  and  O.J.  Marsh,  Appl.  Phys .  Lett.  34, 

^-^C.E.  Jones,  D.E.  Schafer,  W.  Scott,  and  R.J.  Hager,  Bull.  Am. 
Phys.  Soc.  24,  276  (1979)  BL14 . 

^-3M.C.  Ohmer  and  J.E.  Lang,  Appl  Phys.  Lett.  34,  750  (1979). 
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-v.  Swaminathan,  J.E.  Lang,  P.N.  Heminger,  and  S.R.  Smith,  Appl . 
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7G.W.  Ludwig  and  H.H.  Woodbury  in  Solid  State  Physics,  13,  F. 
Switz  and  D.  Turnbull,  Eds.  (Academic  Press,  New  York,  1962), 
p.  223. 
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J.  Schneider  talk  presented  at  the  International  Conference  on 
the  Physics  of  Pure  Semiconductors,  Stockholm,  Sweden,  3-5 
Sept.  1979. 
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EXPERIMENTAL 


Several  experimental  techniques  were  used  to  characterize  the 
samples  and  to  determine  the  properties  of  the  X-center.  These 
include  low-temperature  infrared  absorption,  thermal  annealing, 
deep-level  transient  spectroscopy  (DLTS)  and  Hall  analysis. 

Infrared  Sepctroscopy 

Transmittance  measurements  were  made  using  a  Digilab  FTS-14  Four¬ 
ier  transform  spectrometer  at  a  resolution  of  2  cm-1.  The  samples 
were  lapped  and  polished  with  about  0.5°  taper  to  eliminate  inter¬ 
ference  fringing  and  cooled  to  about  8K  for  the  measurements. 

The  transmittance  was  converted  to  absorption  coefficient  in  the 
D-9 

usual  way  and  digitally  stored  in  the  spectrometer.  The 
transmittance  of  a  high-purity  vacuum  float-zone-grown  crystal 
measured  and  used  as  a  reference  to  subtract  out  any  lattice 
absorption,  particularly  in  the  400-  to  1000-cm-1  region.  Signal 
averaging  of  1200  to  2400  scans  was  used  to  increase  the  signal- 
to-noise  ratio  to  obtain  sufficient  accuracy  with  the  low  con¬ 
centrations  . 

The  samples  were  taken  from  a  variety  of  sources  and  grown  by  a 
variety  of  different  growth  techniques.  Most  of  the  indium-doped 
crystals  were  grown  by  a  solution-growth  process  developed  in 
our  laboratory ,D  with  a  few  grown  by  the  Czochralski  method, 


T.S.  Moss,  Optical  Properites  of  Semiconductors  (Academic 
Press,  New  York,  1959)  p.  14. 

D-10W.  Scott  and  R.J.  Hager,  J.  Electron.  Mater.  8,  581  (1979). 
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and  one  float-zone-grown  crystal.  The  aluminum-doped  and 

gallium-doped  samples  were  obtained  from  crystals  grown  in  our 
laboratory  by  the  Czochralski  technique.  The  boron-doped  crystals 
were  purchased  from  Dow-Corning  and  consisted  of  crystals  which 
were  grown  by  the  Czochraliski  method  and  crystals  which  were 
grown  by  the  float-zone  method.  The  results  we  obtained  did  not 
appear  to  be  unique  to  any  one  particular  growth  method  nor  to 
any  source  of  crystals.  The  specifics  on  crystal  identifications 
are  not  emphasized  here  for  this  reason,  but  where  appropriate, 
any  further  crystal  identification  is  included  with  the  results 
being  discussed. 


The  concentrations  of  the  group  IIIA  acceptor  in  the  crystals  are 
determined  from  the  room  temperature  Hall  coefficient  for  the 
case  of  indium-doped  silicon,  and  from  the  resistivity  curves  re¬ 
cently  published  by  Sclar  for  boron-,  aluminum-,  and  gallium-doped 

. , .  D-12 

silicon . 


Low-temperature  optical  absorption  shows  a  series  of  peaks  due  to 

transitions  from  an  s-like  ground  state  to  p-like  effective-mass 
D-13 

excited  states.  The  concentrations  determined  from  the  Hall 

and  resistivity  data  hve  been  used  to  determine  a  calibration  con¬ 
stant  relating  the  area  of  the  strongest  of  these  peaks  at  8K  to 
the  concentration  of  the  group  IIIA  element.  The  8K  optical  spectrum 
was  then  used  on  most  samples  to  determine  acceptor  concentration. 


D  11Float-zone  indium-doped  silicon  supplied  by  Dr.  Sumner  of  U.S. 
Army  NVL. 

D-12N.  Sclar,  IEEE  Trans.  Electron.  Devices  Ed-24,  709  (1977). 
D-^2D.  Schechter,  J.  Phys.  Chem.  Solids  23 ,  237  (1962). 
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Both  interstitial  oxygen  and  substitutional  carbon  have  vibrational 
bands  which  can  be  seen  in  the  IR  absorption.  The  relationship  be¬ 
tween  the  peak  absorption  and  the  oxygen  or  carbon  concentrations 
has  been  estimated  at  room  temperature  and  is  given  in  ASTM  F120- 
75  and  F121-76  for  oxygen  and  in  Newman  and  Wills0-'1'4  for  carbon. 
Samples  were  measured  at  room  temperature  using  these  calibration 
factors  to  determine  the  oxygen  and  carbon  concentrations.  The 
same  samples  were  then  measured  at  77K  and  8K  to  determine  low- 
temperature  calibration  factors.  The  wavelengths  and  the  8K 
calibration  constants  are  given  in  Table  D-l  for  the  acceptors 
and  for  carbon  and  oxygen.  The  concentrations  of  the  X-centers 
were  determined  solely  from  optical  absorption  measurements  as 
will  be  described  in  detail  in  the  next  section. 


TABLE  D-l.  OPTICAL  CALIBRATION  FACTORS  AT  8K 


Impurity 

Line  (cm-4) 

-3 

Concentration  (cm  ) 

Boron 

279 

[b]  =  peak  area/6.6  x  10-44 

Boron 

666 

=  peak/6  x  10-15 

Aluminum 

472 

=  area/2.4  x  lO--*-^ 

Gallium 

%  510 

=  area/  2  x  10-14 

Indium 

1175 

=  area/5  x  10--*-® 

Thallium 

1907 

=  area/8.5  x  10-4^ 

Carbon 

600 

=  a  peak  x  6.7  x  10-*-® 

Oxygen 

1136 

=  a  peak  3.09  x  10-*-® 

B-X 

184 

=  area/9.5  x  10“-*- 4 

Al-x 

373 

=  area/4  x  10--*-4 

Ga-X 

381 

=  area/4  x  10“ ^ 

In-X 

830 

=  area/3.5  x  10--*-® 

Phosphorus 

315 

=  area/9.94  x  10--*-4 

D-14 

R.C.  Newman  and  J.B.  Willis,  J.  Phys.  Chem.  Solids  26,  373 
(1965). 
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Thermal  Annealing 


The  X-center  concentrations  depend  on  the  sample's  last  annealing 
temperature.  Before  the  infrared  spectra  are  taken,  the  samples 
are  placed  in  a  furnace  in  dry  Ng  gas  and  annealed  until  equilibri¬ 
um  X-center  concentrations  are  reached.  The  equilibrium  condition 
was  checked  by  annealing  the  samples  for  a  longer  time  at  the 
same  temperature  to  see  if  the  concentration  changed.  For  float- 
zone  silicon,  the  annealing  times  used  were  ^  12  hours  and  near 
600°C  and  6  hours  near  1000°C.  For  Czochralski  silicon  ^  4 
hours  near  600°C  and  ^  2  hours  near  1000°C  were  used.  These  times 
are  conservative  and  are  not  the  minimum  times  needed  for  equil¬ 
ibrium.  Baron  et  al .  used  0.5  hour  and  1  hour  anneals  at  650°C 
and  850°C,  respectively,  for  Czochralski  crystals  ~  but  they 
found  these  times  to  be  inadequate  for  float-zone  crystals.  We 
found  that  a  6-hour  anneal  is  too  short  for  equilibrium  to  be 
reached  at  600°C  for  float-zone  crystals.  Once  equilibrium  is 
reached,  the  X-level  concentrations  in  float-zone  and  Czochralski 
silicon  follow  the  same  curves  as  is  shown  later  in  Figure  D-8 . 

The  time  to  reach  equilibrium  is  much  slower  in  float-zone  crys¬ 
tals  however. 

Deep-Level  Transient  Spectroscopy 

Deep-level  transient  spectroscopy  (DLTS)^  ^  ^  has  been  used 

to  study  trapping  in  indium  and  thallium  doped  silicon.  This 
technique  measures  the  capacitance  of  a  reverse-biased  diode.  The 


D  15D.V.  Lang,  J.  Appl.  Phys .  45,  3014  and  3023. 

D  16 

G.L.  Miller,  D.V.  Lang,  and  L.C.  Kimmerling,  Am.  Rev.  of  Mater. 
Sci.  7,  (1977). 
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capacitance  depends  on  the  fixed  charge  in  the  depletion  region. 
N+p  diodes  were  made  in  silicon  doped  with  boron  as  well  as 
indium  or  thallium.  The  boron  concentration  is  held  to  about  10 
times  that  of  the  indium  or  thallium.  Short  voltage  pulses  are 
used  to  bring  free  carriers  into  the  depletion  region  where  some 
of  them  are  trapped.  After  the  filling  pulse,  these  trapped 
carriers  are  thermally  emitted,  producing  a  transient  in  the 
observed  capacitance.  The  signal  is  processed  electronically  to 
provide  a  "spectral"  display  with  a  peak  for  each  trap.  The  peak 
heights  can  be  related  to  the  trap  concentration  and  the  peak 
position  gives  the  emission  rate  as  a  function  of  temperature. 
The  change  in  the  emission  rate  with  temperature  is  used  to 
calculate  the  trap  energy. 

Hall  Technique 


D-C  Hall  data  were  also  taken  on  many  of  the  samples  over  the 
temperature  range  of  40K  to  300K.  The  Hall  coefficient,  re-, 
sistivity,  and  mobility  were  measured.  A  computer  program  which 
does  a  multivariable  least  squares  fit  to  these  data  was  used  on 
some  samples  to  fit  concentrations  and  energy  levels. 

RESULTS 

In  indium-doped  silicon,  an  indium-dependent  acceptor  level 

called  the  X-level  has  been  observed  in  Hall  data.^*  ^  This  level 

occurs  at  0.113  eV,  72%  of  the  substitutional  indium's  energy, 

-2-4 

and  in  concentrations  v  10  to  10  times  that  of  the  indium 

in  pulled  crystals.  This  level  was  subsequently  seen  in  the  in- 

D— 2 

f fared  optical  absorption  spectra.  Figure  D-l  shows  an  ab¬ 
sorption  spectra  for  Czochralski-grown  indium-doped  silicon 
showing  the  606-cm-'1'  carbon  band,  the  113-cm-1  oxygen  band,  and 
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the  indium  lines  near  1200  cm 


Figure  D-l.  Optical  absorption  spectrum  of  indium-doped 
silicon  at  8K  showing  the  prominent  indium 
lines  in  addition  to  the  lines  due  to 
oxygen  and  carbon . 

Figure  D-2 ,  a  scale  expansion  of  the  region  lower  in  energy  than 

the  indium  lines,  shows  a  weak  absorption  spectra,  labeled  X, 

which  has  the  same  characteristic  four-line  acceptor  effective- 

D_2  D— 13 

mass  spectra  as  indium  has.  ’  This  type  of  spectra  arises 

from  transitions  from  an  s-like  ground  state  to  p-like  excited 
states.  The  s-state  penetrates  the  acceptor's  core  region  and 
differences  in  local  screening  shift  this  level's  energy  differ¬ 
ent  amounts  in  different  acceptors.  The  p-type  excited  states 
have  a  node  at  the  core  and  they  are  diffuse  enough  so  that  their 
energies  are  the  same  as  those  calculated  by  effective  mass 
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theory  for  all  of  the  substitutional  group  IIIA  acceptors.  Thus, 
the  s-state  energy  shift  changes  the  infrared  absorption  energy 
for  various  acceptors,  while  the  constant  p-like  final  states 
give  the  same  characteristic  pattern  to  all  the  acceptor  infrared 
spectra . 


Figure  D-2 .  Optical  absorption  spectrum  of  indium-doped 
silicon  at  8K  showing  the  indium-X  lines 
for  a  sample  with  4  x  10l?  In/cm3. 

The  energy  for  an  acceptor  can  be  found  by  adding  the  energy  seen 
in  the  optical  absorption  spectra  for  the  s-to-p  transition  to 
the  calculated  effective  mass  energy  of  the  p-state  to  the  con- 
tinum.  When  this  is  done  for  the  lines  labeled  X,  an  energy  of 
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0.112  eV  is  obtained,  which  is  the  same  value  as  measured  by 
the  Hall  technique  for  the  In-X  center. 

The  optical  spectra  for  silicon  doped  with  the  other  group  IIIA 
elements  (B ,  A1 ,  Ga,  and  Tl)  have  now  been  studied.  In  every  case 
except  for  thallium,  a  defect  spectrum  is  observed  at  energies 

slightly  less  than  the  substitutional  acceptor's  energy,  with 

-2  -5 

intensities  down  10  to  10  in  magnitude.  These  spectra  have 
been  designated  as  B-X,  Al-X,  Ga-X,  and  In-X  spectra  and  are 
shown  in  Figure  D-3  along  with  the  substitutional  boron  spectrum 
for  comparison.  It  is  seen  that  in  each  case  the  number  of  lines, 
intensity  ratios,  and  line  splittings  are  similar  to  those  ob¬ 
served  in  the  boron  acceptor  effective-mass  spectrum.  The  energy 
differences  are  listed  in  Table  D-2  showing  the  same  character¬ 
istic  spacing  as  the  boron  acceptor  spectrum.  Having  identified 
the  spectra  as  being  acceptor  effective-mass-like  spectra,  the 
X-center  energy  can  be  obtained  by  adding  the  measured  s-to-p 
optical  absorption  energy  to  the  calculated  effective-mass  p- 
to-continuum  energy.  These  energies  are  given  in  Table  D-3.  The 
ratios  of  the  X  to  acceptor  energies  are  seen  to  run  in  a  smooth 
progression  from  0.84  for  boron  to  0.72  for  indium. 

As  was  described  in  the  previous  experimental  section,  Hall  and 
resistivity  data  were  used  to  determine  the  calibration  factors 
so  that  the  group  IIIA  acceptor  concentrations  could  be  found 
from  the  strong  (no.  2)  s-to-p  optical  cross  sections  at  8K. 

These  data  are  plotted  in  Figure  D-4.  The  optical  cross  sections 
for  the  s-to-p  absorption  varies  smoothly  with  the  acceptor 
energy.  Deeper  acceptors  have  smaller  s-like  wave  functions  and 
smaller  cross  sections.  The  curve  drawn  through  the  data  for  the 
substitutional  acceptors  was  used  to  estimate  the  calibration 
factors  for  the  X-centers.  These  are  also  shown  in  Figure  D-4, 
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Figure  D-3.  Absorption  spectra  of  substitutional  boron 
in  silicon  and  of  the  X-centers  observed 
for  boron,  aluminum,  gallium,  and  indium. 


and  the  corresponding  calibration  factors  from  peak  area  to  con¬ 
centration  are  listed  in  Table  D-l.  The  concentrations  for  all 
of  the  X-centers  are  then  calculated  from  the  area  of  the  strong 

(no.  2)  peak  in  the  absoprtion  spectra  at  8K.  The  X-center  con- 

—2  —5 

centrations  observed  run  from  about  10  to  10  those  of  the 
substitutional  group  IIIA  acceptor  concentrations.  * 
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TABLE  D-2 .  OPTICAL  ENERGIES  AND  SPACINGS  FOR  THE 
X-CENTER  AND  BORON 


Energy 

/— N 

o 

3 

l 

h-1 

S; 

pacings  (cm  ; 

Defect 

Line 

Line 

Line 

Line 

Lines 

Lines 

Lines 

1 

2 

3 

4 

1  to  2 

2  to  3 

2  to  4 

Boron 

245 

279 

310 

323 

34 

31 

44 

B-X 

184 

220 

250 

261 

36 

31 

41 

Al-X 

338 

372 

404 

414 

34 

32 

42 

Ga-X 

347 

381 

423 

34 

42 

In-X 

799 

831 

863 

873 

32 

32 

42 

TABLE  D-3 .  X-CENTER  DATA 


Dopant 

E(X) 

(meV) 

E(ACC) 

Substitutional 

Dopant 

(meV) 

E(X) /E( ACC) 

[X]/ [ACC] 
Pulled  Si 

ea 

(eV) 

B 

37.1 

44.3 

0.84 

IQ"2 

0.03 

A1 

56.3 

68.5 

0.82 

10_2-10-3 

0.28 

Ga 

57.0 

72.7 

0.79 

O 

1 

1 

H 1 

O 

1 

CJl 

0.40 

In 

112.8 

156.0 

0.72 

10  2-io  4 

0.  70 

T1 

180.0 

246.0 

0.73 

10-2 

— 
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Figure  D-4 .  Calibration  factors  for  converting  the 
area  of  peak  no .  2  in  the  acceptor 
absorption  spectra  at  10K  to  concentra¬ 
tion.  The  B,  Al,  Ga,  and  In  points  are 
experimental  data.  The  X-center  factors 
are  estimated  from  these  data. 

The  indium-X  center  was  also  observed  in  deep-level  transient 
D-17 

spectroscopy.  A  spectrum  for  indium-doped  silicon  is  shown 

in  part  (a)  of  Figure  D-5 .  Four  hole  traps  are  seen  in  the  spec¬ 
trum.  The  DLTS  energies  were  measured  at  0.095  eV,  0.126  eV, 
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Figure  D-5 .  DLTS  spectra  for  indium-doped  silicon 
(a)  and  thallium-doped  silicon. 


0.31  eV  and  0.45  eV.  Electric  fields  lower  the  observed  energies 
D— 18 

of  acceptors  and  this  was  found  to  be  the  case  for  the  first 

D-17 

two  traps.  A  correction  for  the  electric  field  effects  de- 

D— 18 

rived  by  Hartke  was  applied  to  these  two  traps.  The  corrected 

energies  are  0.11  eV  and  ^  0.143  eV  near  the  proper  indium-X 
and  indium  energies  of  0.113  eV  and  0.15  eV.  a  DLTS  spectrum  for 
thallium-doped  silicon  is  shown  in  part  (b)  of  Figure  D-3.  The 


D— 17 

C.  Jones  and  G.  Johnson,  IRIS  Specialty  Group  in  Infrared  De¬ 
tectors,  Minneapolis,  MN  (June  1979).  Also  sub.  J.  Appl .  Phys . 
1980. 

D-18J.L.  Hartke,  J.  Appl.  Phys.  39,  4871  (1968). 
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structure  is  similar  to  that  in  the  indium  spectrum,  with  a  major 
peak  for  the  major  dopant,  thallium,  a  weaker  lower  energy  peak 
to  its  left,  and  two  higher-energy  peaks.  Again,  the  DLTS  energy 
for  the  major  dopant,  thallium,  is  low,  being  ^  0.20  eV  instead 
of  0.245  eV,  and  the  lower-energy  peak  occurs  at  0.14  eV.  If 
a  field  correction  is  assumed  for  the  first  two  peaks  of  0.40  eV 
to  bring  the  thallium  energy  up  to  its  proper  value,  the  peak  mark 
ed  Tl-X  would  have  an  energy  of  ^  0.18  eV.  Hall  and  optical  absorp 
tion  data  were  taken  on  adjacent  pieces  of  thallium-doped  silicon 
crystal.  The  optical  absorption  showed  'v  1016  Tl/cm3  but  no 
shallower  X-like  level.  The  thallium  concentrations  and  energies 
determined  from  the  optical  data  were  put  into  the  computer  pro¬ 
gram  fitting  the  Hall  data  and  the  program  was  allowed  to  deter¬ 
mine  the  parameters  for  any  other  acceptors  needed  to  fit  the 
data.  The  computer  fit  the  Hall  data  with  9.3  x  1015  thallium  and 
another  acceptor  at  0.18  eV  and  4  x  1014/cm3  concentration,  con¬ 
firming  the  DLTS  data  on  the  lower-energy  defect. 

The  Tl-X  center  at  'v  0.18  eV  is  much  deeper  than  the  other  X- 

centers,  giving  it  tighter  wave  functions  and  a  smaller  optical 

cross  section.  This  fact  and  the  low  thallium  concentrations 
16  17  3 

(10  -10  /cm  )  in  the  samples  studied  would  make  it  very  hard 

to  observe  a  Tl-X  absorption  spectrum.  As  yet,  we  have  been 
unable  to  distinguish  an  X-type  spectrum  in  the  thallium  infrared 
data. 

The  concentrations  of  the  X-centers  have  been  observed  to  change, 
depending  on  annealing.  Data  for  the  indium-X-center  concentra¬ 
tion  is  presented  in  Figure  D-6 .  Different  concentrations  are 
seen  in  different  samples  but  in  all  of  these  samples  the  X- 
center  concentration  decreases  as  the  temperature  increases  and 
can  be  fit  to  an  exponential  form^-3, 
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Figure  D-6.  Indium-X-center  concentration  as  a 
function  of  annealing  temperature 
(circles  are  Honeywell  data,  tri¬ 
angles  are  Hughes  dataD~3). 


[in-X]  =  [ln-X]Qe  Ea/KT  (D-l) 

with  an  activation  energy  of  0.70  eV.  The  X-centers  do  not 
anneal  out  up  to  the  silicon  melting  point  and  the  annealing  is 
cyclic.  The  concentration  can  be  lowered  by  annealing  at  a  high 
temperature  and  then  brought  back  to  its  original  value  by  an¬ 
nealing  it  at  its  original  annealing  temperature. 
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The  infrared  absorption  spectra  have  been  used  to  monitor  the  X- 
center,  carbon,  oxygen,  boron,  aluminum,  gallium,  indium,  and 
thallium  concentrations.  The  data  for  indium-doped  silicon  are 
shown  in  Figure  D-7.  The  indium-X-center  concentration  is  normal¬ 
ized  to  the  indium  concentration  and  plotted  with  respect  to 
carbon,  oxygen,  indium,  and  boron.  The  aluminum,  gallium,  and 
thallium  concentrations  were  all  too  small  to  be  detected  op¬ 
tically,  placing  the  electrically  active  concentration  of  alumi¬ 
num  at  less  than  2  x  1012/cm3,  of  gallium  at  less  than  2  x  1012/ 
cm  ,  and  of  thallium  at  less  than  2  x  10l5/cm3.  The  X-center 
concentrations  used  in  Figure  D-7  were  all  adjusted  to  600°C 
equilibrium  values  using  Equation 


Figure  D  7.  Normalized  indium— X-center  concentrations 
( [In-XJ / [In]  )  versus  carbon,  oxygen,  boron, 
and  indium  concentrations.  The  X-concentra- 
tions  have  been  adjusted  using  the  data 
from  Equation  (D-l)  to  600OC  equilibrium 
values . 
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The  data  in  Figure  D-7  for  indium-doped  silicon  show  no  In-X- 

center  dependence  on  oxygen ,  boron ,  or  more  than  one  indium  atom 

per  complex.  There  is  a  linear  relationship  seen  between  the  [X]  / 

[in]  ratio  and  the  carbon  concentration.  The  slope  of  the  line 

through  the  data  is  1,  implying  the  X-center  depends  on  only  a 

single  carbon  per  complex.  The  carbon  dependence  shown  in  Figure 

D-7  is  the  same  as  that  reported  by  the  Baron  et  al.  group  for 

D—3 

the  samples  they  measured.  The  carbon  dependence  was  checked 

by  growing  two  crystals  out  of  indiumD-F9  with  low-carbon,  vacuum 

float-zone  silicon  as  the  source  and  seed  material.  One  sample 

was  intent ially  carbon-doped.  The  two  samples  grown  had  the  same 

indium,  oxygen,  and  aluminum  concentrations,  with  the  undoped 

15  3 

sample  having  less  than  5  x  10  C/cm  and  no  measurable  X-centers 

13  3  16  3 

(<5  x  10  X/cm  ),  while  the  carbon-doped  sample  had  9  x  10  C/cm 

and  2  x  1015X/cm3  after  a  600°C  anneal.  This  confirms  the  carbon 

dependence  of  the  In-X  center  seen  in  Figure  D-7. 

The  simplest  reaction  between  carbon  and  an  acceptor  is 

C  +  A  |  X  (D-2 ) 

D-19 

The  equilibrium  constant,  K,  for  this  reaction  is  equal  to 

K  - - — - =  k  eAF/kT  (D-3 ) 

[A]  [C]  ° 

where  the  brackets  denote  concentration  and  AF  is  the  difference 
in  the  Gibbs'  free  energy  for  the  reaction  in  Equation*-9-2 ^ .  The 
data  in  Figure  D-6  has  been  renormalized  to  [x]  /  [in]  [C]  to  fit 


D  19R.A.  Swalin,  Thermodynamics  of  Solids,  (John  Wiley  &  Sons, 
New  York,  1962)  Chapter  14. 
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Equation^  These  data,  along  with  similar  annealing  data  for 

the  boron,  aluminum,  and  gallium-X  centers  are  shown  in  Figure 
D-8.  The  fact  that  the  normalization  of  Equation^-2)  has  col¬ 
lapsed  the  indium-X  data  in  Figure  D-6  to  a  single  line  is  another 
strong  piece  of  evidence  supporting  the  dependence  of  the  X- 
centers  on  a  single  carbon  and  a  single  acceptor  atom. 


The  Gibbs'  free  energy  difference  in  Equation  (D-3)  is  AF  =  AH  - 
TAS,  where  AH  is  the  change  in  enthaply  and  AS  is  the  change  in 
entropy . 

-AS/  AH/kT 

K  =  KQe  K  e  (D-4) 

In  this  model  the  slopes  of  the  lines  in  Figure  D-8  can  be  related 
to  AH,  which  would  be  expected  to  be  close  to  the  defect  binding 
energies  (0.03  +  0.05  eV  for  the  boron-X  center,  0.28  +  0.05  eV 

for  the  aluminum-X,  0.40  +  0.05  eV  for  the  gallium-X,  and 

0.70  +  0.05  eV  for  the  indium-X  center).  The  infinite  temperature 

intercept  is  related  to  the  defect  entropy.  Differences  in  vib¬ 
rational  frequencies  or  configuration  can  change  the  intercept 
value . 


Mott  and  Gurney  use  an  Einstein  model  where  each  atom  has  three 
vibrational  modes  of  a  single  frequency The  entropy  per  mode 
is  given  as 


S  •  T_ 
vib 


hv 
In  kT 


(D-5  ) 


D  20N.F.  Mott  and  R.W. 


Gurney,  Electronic  Process  in  Ionic  Crystals. 
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Figure  D-8 .  Normalized  equilibrium  concentration  of  the 
X-centers  as  a  function  of  temperature.  The 
quantity  [x]  /[a][c]  is  equal  to  the  equilib¬ 
rium  constant  for  the  reaction  As  +  Cs  |  X. 


The  lighter  atoms  in  the  group  IIIA  series  would  be  expected  to 
have  higher  frequency  modes  and  higher  infinite  temperature  inter¬ 
cepts.  Boron  is  the  only  acceptor  with  a  radius  smaller  than 
silicon.  Because  of  this,  it  could  sit  off-center  from  the  sub- 
titutional  site,  giving  rise  to  a  configurational  entropy  term 
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and  further  raising  the  boron  intercept.  These  effects  are  seen 
as  a  general  trend  in  the  infinite  temperature  intercepts. 

D-21 

The  tetrahedral  radii  for  the  atoms  in  question  are: 


Carbon 

0.774  8 

Boron 

0.853 

Silicon 

1.173 

Gallium 

1.225 

Aluminum 

1.23 

Indium 

1.405 

Thallium 

1.47 

Carbon  and  boron  are  smaller  than  silicon ,  while  the  other  group 

D-3 

IIIA  elements  are  larger.  Baron  et  al  suggested  that  the  in- 
dium-X-center  binding  energy  could  be  due  to  strain  relief  of  the 
substitutional  indium  by  a  nearest-neighbor  carbon.  The  small 
boron  binding  energy  and  the  increase  in  binding  through  the 
group  fits  this  suggestion.  The  fact  that  the  boron-X  center  is 
bound  at  all  and  the  0.12  eV  difference  in  binding  for  the  alu¬ 
minum  and  gallium-X  centers  with  the  same  acceptor  radii  say  that 
there  is  also  a  small  component  to  the  X-center  binding  energy 
due  to  the  defects'  electronic  bonding. 

Stress  Symmetry 

The  absorption  spectrum  of  the  aluminum-X  center  has  been  studied 
as  a  function  of  uniaxial  stress.  A  full  analysis  of  the  stress 
data  is  being  published  in  an  associated  article.  A  summary  of  the 
stress  data  is  shown  in  Figure  D-9  and  a  review  of  the  restuls  is 
given  here . 


D-21 


J.C.  Phillips,  Bonds  and  Bands  in  Semiconductors , 
Press,  New  York,  1973  p.  22 


(Academic 
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Figure  D-9.  Absorption  spectra  showing  substitutional 
boron  and  the  aluminum-X  center  at  8K 
under  unstressed  and  [lOO],  [ill],  and  the 
two  inequivalent  [110]  uniaxial  stress 
conditions . 


Multiple  lines  are  introduced  into  the  optical  data  by  the  stress 
due  to  two  effects.  The  first  is  the  splitting  of  degenerate 
levels  by  the  stress  and  the  second  is  that  a  complex  defect  can 
have  several  possible  orientations  relative  to  an  applied  stress. 
The  second  effect  is  called  a  configurational  effect.  The  data 
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shown  in  Figure  D-9  can  be  explained  by  a  model  involving  the 
following  assumptions: 

1.  The  X-center  ground  state  is  an  orbital  singlet  with 
spin  degeneracy  g  =  2 . 

2.  The  excited  states  will  split  in  the  same  way  as  the 
excited  states  of  the  substitutional  acceptors. 

3.  The  structure  not  explained  by  assumptions  1  and  2 
is  due  to  configurational  effects. 

Assumption  1  is  due  to  the  fact  that  no  optical  absorption  lines 
could  be  identified  as  being  due  to  a  split-off  ground  state.  A 
substitutional  acceptor  has  a  ground— state  degeneracy  of  four 
but  the  crystal  field  from  any  nontetrahedral  complex  would  be 
expected  to  split  these  into  two  spin-doublet  levels.  Spectra 
taken  up  to  50K  showed  no  sign  of  populating  this  second  level, 
so  the  local  crystal  field  is  presumed  to  have  split  these  levels 
into  two  spin  doublets  more  than  10  meV  apart.  Assumption  2  rests 
on  the  argument  that  the  excited  states  are  very  diffuse  and  have 
a  node  at  the  defect.  Figure  D-3  shows  optical  transitions  from 
a  ground  state  to  the  excited  states.  The  X-center  crystal  field 
has  not  perturbed  the  excited  states,  as  is  shown  by  the  effective- 
mass-like  spectra.  The  stress  effects  for  the  excited  X-center 
levels  should  then  be  the  same  as  for. a  substitutional  acceptor 
like  boron.  The  sample  shown  in  Figure  D-9  has  both  aluminum-X 
and  substitutional  boron  in  it  which  can  be  seen  in  the  same 
region  of  the  spectrum.  The  [100]  stress  data  for  the  aluminum-X 
center  have  approximately  the  same  form  as  that  of  boron.  The 
[ill]  and  [110]  stress  spectra  are  more  complex  than  the  boron 
spectra.  The  extra  lines  in  the  latter  two  stress  directions 
can  be  explained  by  two  inequivalent  orientations  of  the  X-com- 
plex  relative  to  the  stress  direction.  The  only  defect  symmetry 
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in  which  the  possible  orientations  would  make  equal  angles  with  a 
•  [100]  stress  but  two  different  angles  with  respect  to  a  [ill]  or  a 

[llO]  stress  is  a  trigonal  [ill]  oriented  defect.  Each  of  the  four 
equivalent  [ill]  directions  of  such  a  defect  structure  makes  the 
same  angle  with  respect  to  a [100]  stress  but  there  are  two  angles 
formed  between  these  directions  and  a  [ill]  or  [110]  stress . 

DISCUSSION 

A  Class  of  Defects 


In  the  first  part  of  this  paper  a  group  of  defects  in  silicon 
was  singled  out  for  study.  These  defects  occurred  in  silicon 
doped  with  a  group  IIIA  acceptor,  with  one  defect  being  associat¬ 
ed  with  boron,  aluminum,  gallium,  indium,  and  thallium.  The  con- 

-2 

centrations  observed  in  Czochralski-grown  silicon  are  ^  10  to 

_5 

10  those  of  the  corresponding  group  IIIA  element's  concentra¬ 
tion.  The  optical,  Hall,  and  DLTS  data  have  shown  that  these 
centers  are  electrically  active  acceptors  having  energies  ^  75% 
to  85%  that  of  the  respective  substutional  group  IIIA  element's 
energy.  All  of  these  centers  showed  the  same  high-temperature 
stability  and  cyclic  annealing  properties.  Taken  as  a  whole,  the 
similarities  in  properties  for  the  defects  we  have  called  X- 
centers  strongly  identifies  them  as  being  a  class  of  similar  de¬ 
fect  complexes.  Specifically: 

•  Each  X-center  is  seen  in  silicon  doped  with  a  different 
group  IIIA  acceptor. 

•  The  X-center  energies  are  all  'v  80%  those  of  the  cor¬ 
responding  acceptor's  energy. 
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•  They  are  all  stable  up  to  the  silicon  melting  point, 
with  larger  concentrations  resulting  from  lower- 
temperature  annealing. 

•  They  can  be  cyclically  annealed. 

•  The  annealing  kinetics  all  fit  A  +  C  +  X  reaction 
kinetics . 

•  The  binding  energies  follow  a  smooth  progression 
through  the  group. 

X-Center  Model 


Aluminum  has  been  suggested  as  being  involved  in  the  In-X  cen- 
D-5 

ter.  If  this  is  true,  we  feel  substitutional  aluminum  should 

have  been  seen  after  high-temperature  anneals  when  the  X-center 

15  3 

is  broken  apart.  The  concentrations  expected  are  in  the  10  /cm 
range  for  some  samples  and  this  is  far  above  the  possible  com¬ 
pensation  levels  which  could  make  it  go  undetected.  Iron  has 

D— 8 

also  been  suggested.  So  far,  we  have  not  seen  iron  in  the 
DLTS  data  but  iron  cannot  be  eliminated  on  the  basis  of  our  data. 
The  chance  exists  that  iron  is  completely  tied  up  in  complexes 
that  we  cannot  detect . 

Figure  D-7  shows  a  linear  relationship  between  the  indium-X- 
center  concentration  and  carbon.  No  relationship  is  seen  for 
boron,  aluminum,  or  oxygen.  Figure  D-8  shows  that  normalizing  the 
indium-X  concentration  to  both  indium  and  carbon  collapsed  the 
indium  annealing  curves  in  Figure  D-6  to  a  single  curve.  This 
normalization  fits  the  form  predicted  from  the  thermodynamics  of 
the  reaction  A  +  C  |  X,  where  K  =  KQexp(EA/KT)  =  [x]  /  [A][C]  .  Two 
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indium-doped  crystals  grown  from  the  same  materials,  except  that 
one  was  carbon-doped,  had  no  measurable  In-X  in  the  undoped  cry- 
stal  and  2  x  10  in-X/cm  in  the  carbon-doped  crystal.  All  of 
these  factors  identify  the  X-centers  as  involving  carbon  and  a 
group  IIIA  acceptor. 

The  irradiation  data u~  suggested  that  vacancies  or  interstitials 
may  be  involved  in  the  X-centers.  The  annealing  temperatures  for 
some  of  the  known  acceptor  interstitial  or  vacancy-related  de¬ 
fects  are  given  in  Table  D-4.  These  annealing  temperatures  are 
all  within  a  few  hundred  degrees  of  room  temperature.  In  none 
of  these  cases  can  the  defect  concentrations  be  cyclically  an¬ 
nealed.  Once  a  vacancy  or  interstitial  center  is  annealed  to  a 
low  concentration,  reannealing  at  a  lower  temperature  does  not 
bring  back  the  original  concentration.  It  is  the  high  thermal 
stability  and  cyclic  annealability  of  the  X-centers  that  leads 
us  to  the  conclusion  that  they  involve  substitutional  atoms  and 
not  interstitials  or  vacancies. 

The  X-centers  have  a  binding  energy  tending  to  form  the  complex, 
while  thermal  energy  tends  to  break  the  components  apart .  This 
fits  the  concentration  versus  annealing  data  shown  in  Figure  D-8. 
Again,  if  the  components  of  the  complex  are  substitutional  atoms, 
they  are  not  eliminated  when  the  complex  is  thermally  broken  a- 
part .  This  explains  the  ability  to  cyclically  break  up  and  then 
reform  the  complexes  when  the  sample  is  brought  to  equilibrium 
at  high  and  then  low  temperatures. 

The  X-center  model  suggested  by  the  data  now  involves  a  substi¬ 
tutional  acceptor  and  a  substitutional  carbon  atom.  It  is  the 
stress  data  which  puts  these  atoms  at  nearest-neighbor  sites.  The 
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TABLE  D-4 .  ANNEALING  TEMPERATURES  FOR  SOME  INTERSTITIAL 
AND  VACANCY-RELATED  DEFECTS 


Defect 

Approx.  Annealing  Temperatures3- 
(time  =  1000  sec) 

V° 

150K 

V= 

7  OK 

P-V 

400K 

As-V 

440K 

Sb-V 

460K 

Al-V 

500K 

A1  -Al. 

S  1 

520K 

C  -C. 

S  1 

550K 

Al. 

l 

500K 

V-V 

550K 

o-v 

600K 

Si 

Data  complied  by  H.J.  Stein,  Radiation  Effects  in  Semi¬ 
conductors  (Eds.  J.W.  Corbett  and  G.D.  Watkins,  Gordon 
and  Breach,  London  1971)  pp.  125-140,  and  R.C.  Newman, 
Infrared  Studies  of  Crystal  Defects  (Taylor  and  Francis, 
London,  1973). 

nearest-neighbor  axis  in  silicon  is  in  the  [ill]  directions ,  which 
was  the  symmetry  found  from  the  stress  data.  The  model  is  shown 
in  Figure  D-10.  Distant  pairs  would  not  be  expected  to  show  the 
distincy  [ill]  stress  symmetry. 
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[too] 


Figure  D-10.  Model  of  the  indium-X  center  showing  the 
stress  axis.  The  carbon  atom  can  occupy 
nearest  neighbor  sites  1,  2,  3,  or  4, 
giving  rise  to  the  different  configura¬ 
tions  under  stress. 


X-Center  Energy  Level  Calculation 

To  determine  if  the  nearest -neighbor  substitutional  acceptor-carbon 
pair  is  a  reasonable  model  for  the  X-centers,  a  calculation  was 
made  to  estimate  this  model's  expected  energy  levels. 

A  mathematical  technique  which  has  been  used  to  obtain  reasonable 
results  for  the  energy  levels,  wave  functions,  and  configurations 
for  defects  in  diamond  and  silicon  is  the  extended  Hiickel  cluster 
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D_22  _j)_20 

technique.  A  defect  is  surrounded  by  three  shells  of 

silicon  atoms  in  the  normal  crystal  lattice  positions.  The  cluster 
is  then  treated  as  a  large  molecule  and  a  liner  combination  of 
atomic  orbitals  (LCAO)  technique  is  used  to  calculate  the  energy 
levels  and  wave  functions. 

The  molecular  orbitals,  Y± ,  are  taken  as  a  linear  combination  of 
Slater  functions,  <f>v,  representing  the  valence  atomic  orbitals  of 
the  cluster  atoms : 

(D-6 ) 
(D-7) 

where  is  a  normalization  constant,  £  is  an  orbital  exponent 
related  to  the  atomic  screening,  n  is  the  principal  quantum  num¬ 
ber,  and  Yln  is  a  spherical  harmonic.  The  energies  of  the  mole¬ 
cular  orbitals  are  calculated  as 

<Yi|H|y±> 

E.  =  - 

1  <Y. I Y.  > 

1  '  1 


D-22 

R.P.  Messmer  and  G.D.  Watkins,  Phys.  Rev.  B7 ,  2568  (1973). 

D-23 

C.  Weigel,  D.  Peak,  J.W.  Corbett,  G.D.  Watkins,  and  R.P. 
Messmer,  Phys.  St at.  Sol,  b  63,  131  (1974). 

D-24 

V.A.  Singh,  C.  Weigel,  J.W.  Corbett,  and  L.M.  Roth,  Phys. 
Stat.  Sol,  b  81 ,  637  (1977).  — 

D-25 

K.L.  Yip,  Phys.  Stat,  Sol,  b  66 ,  619  (1974). 

D  26C.E.  Jones,  Bull.  Am.  Phys.  Soc.  22,  382  (1977)  E014 . 


Y .  =  E  C  .  <b 
1  v  viTv 

-  N.  rn-1e 


In 
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where  H  is  the  Hamiltonian  operator.  The  energies  are  minimized 
with  respect  to  the  coefficients,  C  giving  rise  to  the  secular 
equations 


Z 

v 


E.S  ) 
1  yv' 


C  .  =  0 
vi 


where 


(D-9 ) 


(D-10) 


and 


<t>  > 
v 


(D-ll ) 


The  extended  Huckel  method  involves  approximating  the  H^v  matrix 
elements  by 


(D-12 ) 


K  (H 


yy 


+  H  ) 

vv 


yv 


(D-13 ) 


where  1^  is  an  emperical  ionization  energy  and  K  is  a  constant 
taken  to  be  1.75.  The  secular  equations  (D-9)  are  solved  for  the 
molecular  orbital  energies,  E^,  and  the  molecular  orbital  coef¬ 
ficients,  Cv^ .  The  total  energy  in  the  extended  Huckel  technique 
is  taken  as  being 


E ,  ,  ,  =  Z  n  .  E . 

total  .  l  l 

l 


where  n^  is  the  occupation  number  of  the  ith  state. 


(D-14 ) 
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The  calculations  were  first  run  on  some  known  test  cases:  pure 

silicon,  substitutional  indium,  substitutional  carbon,  and  inter- 

stital  carbon.  The  parameters  used  are  given  in  Table  D-5.  The 

results  are  shown  in  Figure  D-ll.  The  pure  silicon  calculation 

gave  a  bandgap  of  1.01  eV,  with  E  7.4  eV  below  the  vacuum  level 

c 

and  with  an  11-eV-wide  valance  band.  This  compares  well  with  the 

experimental  values  of  1.16  eV  for  the  bandgap,  with  E  5.4  eV 

c 

below  the  vacuum  level  and  a  13.6-eV-wide  valance  band.  Carbon 
substitutional  gave  no  defect  levels  in  the  gap  as  is  the  ex¬ 
perimentally  observed  case.  Carbon  interstitial  in  the  split  [100] 

D-27 

configuration  suggested  by  EPR  results  gave  a  deep  donor 

level  as  is  seen  experimentally.  Substitutional  indium  gave  a 

calculated  level  at  E  +  0.166  eV  close  to  the  actual  level  at 

v 

E  +  0.156  eV. 
v 

The  position  of  a  defect's  energy  level  has  been  found  to  depend 

on  the  interaction  of  the  defect's  atomic  levels  and  the  levels 

D— 28 

of  the  same  symmetry  in  the  host  crystal's  bands.  The  good 

results  for  the  pure-silicon  bandgap  and  bandwidths  and  the  good 
results  for  the  indium  and  carbon  test  defects  suggest  that  as 
crude  as  the  extended  Hiickel  approximation  is,  it  may  be  cal¬ 
culating  a  reasonable  average  density  of  states  for  the  silicon 
crystal.  This  gives  some  reason  to  suggest  that  the  technique 
should  be  able  to  be  used  to  describe  other  similar  defects. 


D— 97 

G. D.  Watkins  and  K.L.  Brower,  Phys .  Rev.  Lett.  36,  1329  (1976). 

D— 2  8 

H.  Hjalmarson,  Thesis,  U.  of  Illinois,  Physics  Dept.  (1979). 
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TABLE  D-5 .  PARAMETERS  USED  IN  EXTENDED  HUCKEL 

CALCULATIONS 


Atom 

Wave 

Function 

Energy 

(eV) 

Exponent 

Silicon 

3s 

-14.66f 

1 . 6344a 

3p 

-  8.068 

1 . 4284a 

Surface  silicon 

3s 

-  8.1 

1 . 63a 

Indium 

5s 

-11 . 47C 

1 . 902 3b 

5p 

-  5.787 

1 . 6940b 

Carbon 

2s 

-19 . 20e 

1 . 6083a 

2p 

-11.79 

1 . 5679a 

Aluminum 

2s 

-10 . 62d 

1 . 3724a 

2p 

-  5.986 

1 . 3552a 

Oxygen  ^ 

2s 

-13.62 

j  2.2458 

2p 

-17.19 

2.22 66a 

aE.  dementi  and  D.L.  Raimondi,  J.  Chem.  Phys. ,  38, 
2686  (1963). 

°E.  Clement i,  D.L.  Raimondi,  and  W.P.  Reinhardt,  J. 
Chem.  Phys.  47,  1300  (1967). 

cC.  Moore,  Atomic  Energy  Level  VIII  Circular  of  the 
NBS  No.  467  (1958). 

dW.  Lotz ,  J.  Opt.  Soc.  Am.  60,  206  (1970). 

a  eS.  Huzinaga,  "Approximate  Atomic  Functions." 

f 

G.D.  Watkins,  private  communications. 
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Figure  D-ll.  Energy  levels  calculated  for  pure  silicon, 
substitutional  carbon,  interstitial  carbon, 
and  substitutional  indium  to  test  the  ex¬ 
tended  Huckel  cluster  technique  shown  in 
the  top  diagrams.  The  bottom  diagrams  are 
for  the  indium-X-center  model  of  a  nearest- 
neighbor  substitutional  carbon-indium  pair 
in  the  substitutional  sites  and  in  a  re¬ 
laxed  minimum— energy  configuration. 

The  extended  Huckel  theory  does  not  include  charge  interations, 
and  clusters  with  atoms  having  differences  in  Pauling  electroneg¬ 
ativity  greater  than  1.2  have  been  found  to  give  poor  results.0-29 
Silicon,  indium,  and  carbon  have  electronegativities  of  1.8,  1.7, 
and  2.5,  all  within  0.8  of  each  other. 


D-29 


R.T.  Sanderson  in  Chemical  Periodicity  (Reinhold  Publishing 
Co.,  New  York,  1969)  p.  34. 


4 


With  this  background,  the  nearest-neighbor  substitutional  carbon- 
acceptor  pair  model  for  the  indium-X  center  was  then  calculated. 

A  single  energy  level  slightly  shallower  than  that  of  the  sub¬ 
stitutional  indium  was  found.  The  total  energy  of  the  complex 
can  be  reduced  by  letting  the  atoms  relax,  with  the  indium  and 
the  carbon  moving  away  from  each  other.  In  this  configuration 
the  energy  level  for  the  defect  is  0.13  eV,  close  to  the  actual 
X-center  energy  of  0.113  eV. 

The  wave  function  corresponding  to  the  0.13  eV  level  is  a  bonding 
molecular  orbit  between  the  indium  and  the  three  neighboring 
silicons.  The  carbon  neighbor  does  not  take  part  in  this  bond. 
This  explains  why  the  complex  relaxes  to  a  minimum-energy  con¬ 
figuration  by  moving  the  indium  and  carbon  away  from  each  other 
and  closer  to  their  bonding  neighbors.  The  energy  levels  for  the 
carbon  acceptor  pair  are  also  shown  in  Figure  D-ll. 

The  fact  that  the  extended  Hiickel  cluster  calculation  came  out 
with  an  energy  level  for  the  carbon-indium  pair  which  is  close 
to  the  measured  In-X-center  energy  gives  some  further  support 
to  the  model.  In  this  case  it  is  probably  of  more  interest  to 
take  a  different  view.  That  is  that  the  experimental  data  of 
the  dependence  on  carbon  and  an  acceptor,  the  high-temperature 
cyclic  annealing,  and  the  [ill]  stress  symmetry  have  established 
the  X-center  model  and  that  the  calculations  give  us  some  in¬ 
sight  into  this  defect's  bonding,  wave  functions,  and  minimum- 
energy  configuration.  The  combination  of  the  experimental  data 
and  the  theoretical  calculation  gives  a  good  total  picture  of 
the  carbon-acceptor  pair  (or  X-center)  defects. 
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Application  Note 


One  of  the  reasons  for  interest  in  X-centers  was  their  adverse 
effect  on  extrinsic  silicon  infrared  detectors.  With  the  identi¬ 
fication  of  the  carbon  dependence  for  these  centers,  it  has  been 
possible  to  produce  low-carbon  detectors  with  no  measurable  X- 
level  concentrations. 

CONCLUSIONS 


The  data  for  the  X-centers  in  silicon  yields  two  major  results: 

1.  A  class  of  defects  called  X-centers  exists.  There  is 
one  of  these  centers  corresponding  to  each  of  the 
group  IIIA  elements  boron,  aluminum,  gallium,  indium, 
or  thallium  in  silicon. 

2.  The  model  for  these  centers  which  best  fits  the  data 
is  a  nearest-neighbor  substitutional  carbon-acceptor 
pair . 
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